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A CORRESPONDENCE BETWEEN NULL AND 7 
ORDINARY GEODESICS 


By P. G. GORMLEY (Dublin) 
[Received 23 October 1933] 


THE properties of spaces with corresponding geodesics are well known. 
When a certain coordinate-system is chosen in one space, then there 
exists a coordinate-system in the second space such that the Chris- 
toffel symbols are related by the equation 

feel ec (1) 
where y,, is a gradient vector. 

In this paper it is proposed to determine the relation existing 
between two spaces V, and ‘V, when the null geodesics of V,, corre- 
spond to a set of ordinary geodesics of ‘V,. This is expressed analyti- 
cally as follows: in a certain coordinate-system in V,, the equations 
of the null geodesics are . 

zt = F(t; a,, dg,..., gn), (2) 
where the parameter ¢ alone varies along a particular null geodesic, 
and the 2n—1 parameters 4a, d,... determine this curve uniquely. 
Then in some suitable coordinate-system the equations (2) define 
a system of oo?”-% geodesics in ‘V, and we propose to determine the 
relation existing between the spaces V, and ‘V,. 

1. In general, however, it is not possible to determine the null 
geodesics as parametric equations of the form (2), and then we must 
consider the differential equations which determine them. When the 
metric of the space V,, is 

ds* = g,,, dxtdx”, (1.1) 
then the differential equations* of the null geodesics are, dA being 
the element of spatial distance, 

ial (1) = = , (1.2) 
dd? * \xB} dA dr AC 
da dxB _ SS 
IoB Ty aA nl 


* The definitions of the quantities (a) Ryo, ete., used in this work are 


those given by Eisenhart, Riemannian Geometry. 
3695.5 B 





2 P. G. GORMLEY 
When we use a general parameter ¢ other than A, these equations 
become 

Gl ae , { e\ dam =) dat (in + fv \ dx 


dt \ dt T |aB} dt dt} = dt \ dt? ' (BJ dt dt 
da® dxB _ 


= 0. iB 
dt dt ( 


In a ‘V, with metric d's? = "9, datde’, (1.4) 


=~ QS 


Iup 


the equations of the geodesics are similar to (1.3), except that the 


Christoffel symbols { #\ must be replaced by ‘{#\ 

ap} “ lap) 
formed from the fundamental tensor ’g,,,. Then, if every solution of 
the equations (1.3) and (1.3’) is to define a geodesic in ‘V,, we see 
that the equations 


(1a tallS-[o)- lle) ¢ 2 


must be satisfied when the equation (1.3’) is satisfied. Since 


, these being 


‘(H\_ fe) —_ 7" + 
lap) lap ~~ °F” 

where 7',2" is a tensor, the equations 

la® dx? dav 

5” TT: p—.§4 J: ov ¢ pecnt a. <2 

Oy Sap Oy of) ay at ai 


must be satisfied in consequence of (1.3’). Introducing a Lagrangian 


ee : da 
multiplier, which must be of the form A,” a? we must have 
cf 


identically 
‘ la® dxP day 
8” Tt — Se Tz” —A-"¥g,2)—— —- —- = 0, 1.5 
( Y xB Y ap Y Iup) dt dt dt ( >) 
where A,” is skew in vy and yw. From (1.5) we find 
oY tia — 8, T xh +84 74 —% Tp +8 Px — 85 Sil + 
¥ +AygagtAsgp,+Agg,, = 90. (1.6) 
Summing with respect to v and y we find 
_ 1 1 Quy | A ons ae iad 
7 xp —2 +1 (dé Th+5% T+ A! Jap tA pi +Ag,), (1.7) 
where 7; = T'3,*, and A¥ = A,. 
Then, if we substitute from (1.7) in (1.6), we find an equation 
involving A;”” only. Multiplying by g®v and summing for f and y, 





ON NULL AND ORDINARY GEODESICS 
it is found that, when the suffixes p and v are lowered, 
[(m +4)(n see 1)+2JA guy t 2A pygt2A yay t(n+4)(A, Iva—AyIpa) = 0. 

(1.8) 
Writing the two similar equations found by permuting cyclically apr, 
and adding the three equations thus found, we see that 
(n+-1)(n+-2)(Aguy+Ayvat+Arvay) = 0. 

From this and (1.8) we find 


1 
Bug a n—1 (A, 9ap—A, Joy)s 


which, when substituted in (1.7), gives 
T ph = KebptBebatGapoe, | (1.9) 
1 1 1 
yher PY om Av = —_(T,— -A,|, 1.10 
where dp POS pa Po =al . ae .) (1.10) 
and ~ AzHY = 5h hp —8y gr. (1.11) 
If we substitute the values of 7',g4 and A,” found from (1.9) and 
(1.11) in (1.6), this equation is satisfied identically. Thus the con- 
ditions ' 
he hs 
\wAl (pa 
are necessary and sufficient that the null geodesics of the space V,, corre- 
spond to a set of geodesics of the space 'V,. : 
If the geodesics of ‘V,, are null geodesics, then \ : , 
] 2 
dx®* db G4 
y ~ —— == © 1.13 
IB at di ne 
must be a consequence of (1.3) and (1.3’). So since (1.13) does not 
involve d*x%/dt?, then 


+8n Prt fy t+Iprh” (1.12) 


"Jap = "Gaps “ (1.14) 
i.e. 'V, is conformal* to V,. In this case we have : 

dy = —Yy x = fF, (1.15) 

When the multiplier A,;"” is identically zero, then ¢ is identically 

zero, and conversely. Then we see that all the geodesics of V, corre- 

spond to those of ‘V, and the spaces are in projective correspondence. 

Then 


{A = [teeta te (1.16) 


\uA} pal 


a 9 . 
and $, = 73 log, [2 (1.17) 


is a gradient. 


+y 


* Cf. Eisenhart, Riemannian Geometry, 89. 











4 P. G. GORMLEY 
In the general case we find from (1.12) that 
é 9 
Aclog, [2 = (n = +1)p,.+dy» (1.18) 
and so it is evident that each of the vectors ys, and ¢, need not be 
a gradient. However, it will be seen that the most important cases 
are subject to this last condition. 


2. The element of spatial distance along a null geodesic is a scalar 
multiple of the corresponding element of arc-length of the corre- 
sponding geodesic. It will be seen that this multiple factor depends 
only on the vector y,,. 

The equations of the geodesics of 'V,, are 
@ar '[v\ der det _ 4 
d's” ap) d's d's’ 
d*x” | {v \ dx dab _ 9,,, 1x% da” " dx” dxP 
dst \ap| de ds ds dso 8 as de 
Those geodesics corresponding to null geodesics of V,, satisfy 
dx” {v\ dx dab 9, a% dx” 
di \ap| de d's as de’ 
da* dx _ 
Iop 7 
dx” | {v\ dx dak 





or 


and hence 





dd" \aB) dA dX’ 
’ da* dab _ 
a on 
dr/d's® da 
rher ’ ‘ _s 
neve daA/d's "ba Ts d's” 
Thus = = = exp(—: 2| bm dl : vs), (2.1) 


If the ratio is to be a function of position, then it is necessary and 
sufficient that %,, be a gradient o/ea". Then 
f= = constant x e~2¥, (2.2) 
and since we can multiply A by a constant factor without changing 
the form of the equations of the null geodesics, then 
da 


= = e-*, (2.3) 
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Thus we see the relation to conformal geometry where we are 
given that the ratio of corresponding non-null elements of length is 
ds/d's = e-¥, in which case the ratio of corresponding elements of 
spatial distance is a ds\2 
—_—_ = —2y = i—-]. 

mn * (7) 
Now we have seen that ¢,,+(n+1)%, is a gradient, and hence so 


is ¢,. Consider the V,, with fundamental tensor 


in = oy, | (2.4) 
then La ~ 8 bb. + 9,9" (2.5) 


so (1.12) gives 


(7) (Flog j 
{a} Pa {a} Bead) +84 Hy) (2.6) 


Thus there exists a projective correspondence between ‘V,, and V,. 
When the ratio (2.2) is constant, then %,, = 0 and 


{ial i a tad" (2.7) 


3. Relation between the corresponding curvature tensors, and 
the case where ‘J, is a flat space 
From (1.12) we have* 
Ryo = Royo t8\(Vp to— Vo Py) —Fi(Vo ta— to fra— Gar Pa) + 
+3(V, ta— Py Pr—Ipr Vu) —IprlVoP? +o 9") + 
+9oV ib" +, $"), (3.1) 


whence 
Ry, _ Ry, +-V,, $,—V) $+ (n— 1)(V,, py— Py Pr—GJpr Pa) + 
+V, PratGy Pr—IpnVoh? +ho 9%): (3.2) 
It is easily proved that the tensor ‘R,,, is symmetric when (1.18) is 
taken into consideration. Multiplication of (3.1) by ‘g;4 = ’g,,g% 
and summation with respect to v and A give 
v= GaVy ¥.—VY, $u—"IualVo p°—, P+ 'GJoalVy pr—f, pr) — 
—"GuolVo $*+bo9") + GoalV, $* +h, $"). (3.3) 


* The symbol V, denotes covariant differentiation with respect to the 
metric (1.1). 








P. G. GORMLEY 
When the space ‘V,, is a flat space, then ‘R'),,, = 0 and 
RR’ \yo = 5\(V, $,.—-% to) +97, ona PurtGpa Qa”—IorQ,.”; 


where Pur = Vu Pr— Pp Pr— 29 P Yas 

and = Vu ¢’+do,9'+ 301, Py. 

Then, multiplying (3.4) by g4 and summing, we find 
RY, = Voy —Vpgt Pa” —85 Py t+ (n—1)Q,"; 

while the condition R,, = R,, and equation (1.18) give 


0= Vo ¥,—‘, Ba» 

. ion ae aan ‘ 1 ne — Ay 
or p, = ep/ex¥, and hence ¢, = 6¢/éx". ahs. 

In this case the ratio of the element of spatial distance to the corre- 
sponding element of arc-length (or spatial distance) of the corresponding 
geodesic in the flat space is always a function of position. 

aS ot ee ates Ce Se P oe 
Now the tensors P,, and Q,,, are symmetric since %, and 4, are 
gradients. Since 
— a] > ! ~ 
iin = Iup ee at Jur es Quy (3.5) 
and, by interchange of v with A and p with w, 
— 2 | 
we = Dru Fin oat Jo Qu Juv Qur» 
we have, by addition, 
— | 
oe —s Ivy a Ly Irn digs Jrw | 
here ek 2 
where Lun = 3(Far— @ya)- 

Then, when > 3, the space is conformal to a flat space. 

The geodesics of a flat space can be represented linearly’in terms 
of a parameter, and hence so can the null geodesics of a space related 
to it by the correspondence (1.12). 

Thus the necessary and sufficient condition that the null geodesics of 
a Riemannian space of n (> 3) dimensions can be represented para- 
metrically by an equation of the form 

x’ = att+be (a, b constants), 
is that the space be conformal to a flat space. 

In this case we can choose the null geodesics of the flat space as 
the set of geodesics. Then by (1.15) ¢, = —y%,. In a system of 
cartesian coordinates the metric of the flat space is 

, m av 
d's? = y,,, davdx’, 
where y,, = 0, » AV} Or y,y+1, » = v, while that of the Rieman- 


nian space is ; ; 
i ds? = e-y _, dada”. 
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The differential equations of the geodesics are 


a otf dat da® dx8 
ds? “ds ds +9 


“Waa” 


and for the null geodesics 


d dat 
ay = @. 
an (e- zi 
By choosing a parameter ¢ such that e~?¥ = dA/dt, these equations 


can be written as dat . 
ts 
whence, on integration, at = att+oe, 
where a and b+ are constants. 
The value of A is given by 
t 
A= | e-24 dt, (3.9) 
0 
where we suppose that the value of A is zero when 2° = 6°. From 
the equation y,g a*a® = 0 and (3.8) it follows that 
rx’ —y,,, bebY 
2Y yy ahd” 
Then (3.9) and (3.10) give, apart from a scalar multiplier, the spatial 
distance between two points on the null geodesic. The transport 


vector at a point is defined to be 


t = “uv (3.10) 


dat 


= aA’ 
Then from (3.8) and (3.9) its value is found to be 
nt = e%ar, (3.12) 
if Tn _ Purat Qua (3.13) 
then from (3.5) it follows that 


v= Iv Dew —Ivw Py t+9ru ee a 


1 
Tan = = Te"9ur- (3.14) 


(3.11) 


Now let ¢, = 0. Then Q,,, = 0, and (3.14) reduces to 


1 ‘“ 
P= 7 Pa Oey (3.15) 


and (3.6) shows that the space is one of constant curvature if mn > 2 








8 P. G. GORMLEY 
Hence, if a space satisfying the above conditions is projectively flat, 
then, for n > 2, it is a space of constant curvature.* 


If 4, = 0, then F,, = 0, and 


1 i 
Qu = - Q.°"Gur- (3.16) 


. 7 : = 1 
Thus when x > 3 the space is of constant curvature Ky = —— Q;*. 
n 


Then , 2 44 
Vi Pat Gu dr = —Kogyr (3.17) 
which can be written in the form 
V,, Vie? = —Kye?g,. 
The solution of this equation ist 
e* = cos(K}s), (3.18) 
where s is the length of the geodesic joining the point 2# to a certain 
fixed point in space. 

In the latter case it can be shown that when n = 2, the space is 
one of constant curvature, as has been proved for the previous case. 

If the space ‘V,, is a space of constant curvature, and we are con- 
sidering the general correspondence (1.12), then 

¥ = Ko(8;, ‘Vro—9e Gru) 
As before, it is easily shown that y,, is a gradient and the space JV, 
is conformal to a flat space, if m > 3. 

In particular, if 6” = 0, then V, is a space of constant curvature.§ 
If %,, = 0, it does not seem necessary that V,, be a space of constant 
curvature. 

4. Is the correspondence (1.12) a reversible one? 
Before considering this question we shall state a result which is 
useful for the solution. T'he most general solution of the equation 
0 = Ftp tgrl’ + Gur P+ "Gur "os (4.1) 
where the rank of the matrices (g,,), ('9,), and ("g,,,) 18 m ( > 2), ts one 
of the following: 
(i) $,=0= 9 = '¢" = "9", 
(ii) $.=9, ‘Gur= yr "Tur= 5G, $’+0'd"+b"d” = 0, 
* Weyl, Géttinger Nachrichten (1921), 110. 
+ Cf. Ruse, Quart. J. of Math. (Oxford) 1 (1930), 146. 


t Cf. Schouten, Der Ricci-Kalkul, 131. 
§ Cf. Schouten, loc. cit., 204. 
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(iii) Py = "Our = ‘Gur ¢’ = 0, ‘p’+c"p" = 0, 
and two similar equations, 
> cs ee 


(iv) $.=99 Ga=UGatm’g.», $= ae 


The proof follows at once from the equation found from (4.1) on 
multiplication by g,, where q,"¢” = 0 and q,’¢” = 0, and summing. 

If we have three Riemannian spaces V,, ‘V,, "V,, such that the 
null geodesics of V, correspond to geodesics of ‘V, and "V,, while 
the null geodesics of ‘V,, correspond to geodesics of "V,, then 


= 
J 


| 
| 


\ 
| 


(a 


1) Peat ohh tad, 
BA) 8 a at at 


= pal Tht Yat "8 


Vv 
ur 

Vv 

A 


whence 
0 = 8) 6, +8, I,+9 ure” — ‘”) + ‘Gur "?¢’, 6,, was: $,— ‘but "Pu 


Thus we see that either Pe. 
Jur = yr 


or “?” = 0. 

The correspondence between V, and 'V, is a conformal one or that 
between 'V,, and "V, is a projective one. 

By making "V,, coincide with V,, it is seen that the necessary and 
sufficient condition that the correspondence (1.12) be reversible, is that 
it be a projective or a conformal one (n > 2). 

If the null geodesics of V,, ‘V,, "V,, correspond to geodesics of 'V,,, 
"V,,, V,, respectively, then it can be shown that one of the following 
statements is true: 

(i) the correspondence between each of the pairs V,, 'V,; 'V,, "Vy; 
"V,,, V,, is a conformal one; 

(ii) the correspondence between each of the pairs is a projective one; 

(iii) the correspondence between one of the pairs is a projective one, 
and that between another of the pairs is a conformal one; 

(iv) there exists a relation of the form 

Jur = VGprtb"9,r 
between the metrical tensors. Then it follows that the tensors ‘g,, and 
"Gur have a set of principal directions common to each other in the space 
V,. Also, the vectors ¢”, ‘¢”, and "$” have the same direction. 


n* 
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‘ON THE NUMBER OF REAL ROOTS OF A 
QUINTIC EQUATION 


By T. W. CHAUNDY (Ozford) 


[Received 15 September 1933] 


Outline of the method 
Tue problem of determining the number of real roots of a quintic 


equation 4, + 5a, 4+ 10a, t?-+ 10a, f+ 5a,t-+a, = 0 (1) 


has been discussed by Hermite,* Sylvester,t Salmon,{ Cayley.§ 
Baker.|| The solution, more particularly as given by the last-named, 
who resumes and extends the work of his predecessors, states neces- 
sary and sufficient conditions in terms of invariants of the quintic 
and is theoretically complete, but the technique is advanced and the 
results heavy if employed for practical purposes. 

On the other hand, Sturm’s method (however admirable in a simple 
numerical case), when applied to the general quintic, leads through 
tiresome calculations to inequalities which are not the simplest pos- 
sible; nor even the fewest possible, for Sturm’s method leads to four 
inequalities, whereas the authors cited above and I are alike content 
with three. 

The method which I describe here is based mainly on graphical 
considerations and indeed finds some of its arguments in an appeal 
to the visual geometry of a figure. It is simpler and sufficient to 
consider (1) in the reduced form** 

§+-10A,#+10A,02+5A,1+A; = 0, (2) 
where now A,, As, Ay, A; are seminvariants{f of the original quintic 
(1). Actually I write (2) in the form 

+-10A, t+ 10A,02+-5yt+a = 0 (3) 
and interpret it as a straight line, depending on a variable parameter 

* Hermite, Cambridge and Dublin Math. J. 9 (1854), 172-217 (198-217). 

+ Sylvester, Collected Works, 2 (1908), 371-4, 375, 418-79, 482-3. 

~ Salmon, Modern Higher Algebra (1885), 240. 

§ Cayley, Collected Works, 6 (1893), 147-96. 

Baker, Proc. London Math. Soc. (2) 6 (1908), 122-40. 
** We write a, t+a, = t’ in (1) and then drop accents. 
tt The criteria obtained by this method are stated in terms of seminvariants 


only, like the standard criteria for the quartic, not in terms of full invariants 
like those discussed by the authorities quoted. 
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, in the plane of cartesian coordinates (x, y). Keeping A,, A; con- 
stant, I trace the envelope of this line for variation of the parameter 
and determine the number of real roots of the quintic for given x, y 
as the number of real tangents which can be drawn to the envelope 
from the corresponding point (x, y). 

The envelope itself is the locus obtained by equating to zero the 
discriminant A of (3). As we know from algebraic considerations, 
the inequalities A > 0, A < 0 distinguish the respective cases of an 
even number and an odd number of conjugate pairs of imaginary 
roots. Actually the curve A = 0 does more and completely delimits 
the several regions in the plane in which the quintic has respectively 
0, 2, 4 imaginary roots. The main analytical problem lies in finding 
convenient algebraic definition of these several regions, which have 
the form of curvilinear polygons bounded by finite arcs of the 
envelope. 

The envelope, which is a unicursal curve, can only divide the plane 
into these three or more regions, if it cross itself once or more. The 
determination of the double points of the envelope and in particular 
the discrimination between nodes, cusps, and isolated points becomes 
accordingly important. 


Application to the quartic 
[ shall illustrate the method by applying it in the first instance to 
the general quartic,* obtaining the well-known conditions for 4, 2, 0 
real roots. 


Take the quartic in its reduced form 


f(t) = 446A, ?+4yt+2 = 0. 


Differentiation gives for a typical point of the envelope 


x = 3(t#4+2A,t%) ) 
y = —(#8+34,t)) 
Differentiating along the envelope we have 
dx = 12t(2+.A,) dt, dy = —3(t®?+A,) dt. (6) 


* A similar geometric discussion of the quartic is given in Aubert et Papelier, 
Exercices d’algébre, d’analyse et de trigonométrie (Paris), 2 (1910), 321-3. The 
form of the envelope is determined as in Fig. 2 below, but no attempt is made 
to deduce algebraic criteria from it. Since I regard this as an essential part 
of the method, I have preferred to discuss the quartic de novo. For this 
reference I am indebted to Mr. W. L. Ferrar. 
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The cusps are given by dx = 0 = dy, i.e. by 
f'(t) = @+4A, = 0. (7) 
The curve goes through the origin (at ¢ = 0) and cuts y = 0 again 
where #434, = 0. (8) 
These two values of ¢ give only a single value of x, namely x = 943. 
The point (9A3,0) is therefore a node of the envelope. This gives 
three double points. By (4), (5) the class and degree of the envelope 
are alike four, and so, by Pliicker’s equations, there are no other 
double points. We may add here that there are no finite points of 
inflexion, since, by (4), the ‘slope’ of the 
tangent varies inversely as ¢ and the tan- 
gent therefore continues steadily rotating 
as ¢ moves steadily along the curve. 
We must now distinguish certain cases. 
If A, > 0, the roots of both (7) and (8) 
are imaginary. The cusps are therefore 
imaginary and the node an isolated point. 
The coordinates x, y are both monotonic 
functions of ¢ and the envelope is ‘para- 
bolic’ in shape as in Fig. 1. 


Evidently two tangents can be drawn 
to the curve from points on the convex 
side and none from points on the con- 
cave side. The numbers of real roots are then as marked in the 


Fic. 1. (A, > 0). 


figure. 

If, on the other hand, A, < 0, say A, = —c?, then the cusps are 
at the real points (—3c*, + 2c*) and the nodal tangents are also real. 
The envelope has now the form shown in Fig. 2, and the numbers 
of real roots are again as marked in the figure.* 

There is also the intermediate limiting case in which the cusps, the 
node, and the region [4] of Fig. 2 have moved up to coalescence at 
a single point (the origin) on the curve of Fig. 1. The figure is not 
visually distinguishable from Fig. 1, and need not be specially drawn. 

Each arc of the envelope itself, in either figure, must be ascribed 
to the region on its convex side, i.e. to the region of the greater 


* We can go further and, as in Sturm’s theory, determine the number of 
roots in an interval (a, b) by drawing the tangents at ¢ = a, t = b as further 
boundaries of regions. 
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number of real roots. The points of the envelope correspond to 
repeated roots; the node to two distinct roots and each cusp to a 
triple root. 

If we exclude the limiting case A = 0 for the moment, then, as 
we have already said, the sign of A distinguishes the case of 2 real 
roots from the case of 0 or 4 real roots, and it is in the separation 
of these two latter possibilities that the crux of the problem lies. 
The argument is most easily phrased, if we set out first to obtain the 
exact conditions for all roots real. 


(- 
——K——aee ween ewe 


0 


Fie. 2. (A, < 0). 


Evidently, except for the limiting case A, = 0, A = 0, x = 0, we 
need A, < 0, in order that the envelope may have the form of Fig. 2 
and then, in that figure, A > 0. It remains to separate the regions 
[4] and [0]. Evidently any curve through the node will do this, pro- 
vided that it does not actually cross the envelope again. Let us take 
the simplest such curve, namely, the ordinate through the node, 
shown in Fig. 2 by the broken line. We can then read off from the 
two figures the following exact criteria for real and imaginary roots 
(repeated roots included): 

[4] Four real 

A>0, A, < 0, x < 9A?; (9) 
A= 0, A, = 0, 2 = 043; 
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[2] Two real, two imaginary 
A< 0; 
A = @, A, > 0; 
A, = @, x 4 9A}; 


A,<0, 2x>9A3;) 


[0] Four imaginary 
A ><; A, > 0; ) 


a > 9A?) 


In the standard notation for the quartic* 
A, = H =a,a,—a; 
y = G = apa,—3a,a,a,+ 2a} 
x = a; I—3H?, where I = a,a,—4a,a,+ 3a} 
Thus the conditions (9) give, in the standard notation, 
H< 0, A>0, a; I—12H? < 0, (13). 


which are the standard conditions for four real roots.t (I omit the 
trivial quartic 4 = 0.) 

It is of interest to remark that Sturm’s method gives, in place of 
the third inequality of (13), 


3(G?+ 8H?) > 2a5 HT, 
i.e., in the notation [ am using, 
3y?—2A,x%+18A35 > 0. 


The left-hand side equated to zero gives a parabola through both 
cusps and the node, shown in Fig. 2 by a dotted line. We. could 
evidently have obtained Sturm’s condition by taking this parabola, 
instead of the ordinate x = 9A3, as the curve separating the regions 
[4] and [0]. A perverse ingenuity will suggest even more elaborate 
curves of separation with corresponding inequalities still less advan- 
tageous.{ 


* See, for instance, Burnside and Panton, Theory of Equations, 1 (1904), 73-5. 

+ Burnside and Panton, loc. cit., 144-5. 

t It may be recorded that a similar element of choice in the criteria of the 
quintic is greeted by Sylvester [loc. cit., 463] with characteristic enthusiasm. 
The geometrical approach makes it clear why such elements of choice are to 
be expected in the criteria of reality of roots of any algebraic equation of 
degree exceeding three. 
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The general quintic: cusps and nodes 
I come now to the general quintic in its reduced form (3) 
4-104, 8+10A,0+5yt+2 = 0. 
Differentiation in t gives for the typical point of the envelope 
a = 41+20A,+10A,? ) 
y = —t—6A,t2—4A,t | 
At a cusp dx = 0 = dy, ie. 
20(t4+-3A, 12+ Agt) dt = 0 = —4(®4+3A,t+A,) det. 
Thus the parameters of the cusps are the roots of the cubic 
#+3A,t+A, = 0. (15) 
The discriminant of this ‘reduced’ cubic is known to be —(A3+4A3). 
Thus the cusps are all real, if and only if 
A3+4A?2 < 0. (16) 
[t is of course incidentally necessary, but not sufficient, that A, < 0. 
Now a node corresponds to a pair of double roots of the quintic. 
Thus, if (#,y) is a node at which the two values of the parameter 
are given by #2 2ot+-B = 0, (17) 
we have an identity of the form 
®+10A,8+10A,02+5yt+a = (t+4a)(2—2at+f)?, (18) 
the residual factor on the right being chosen to secure the disappear- 
ance of the term in ¢*. Equating coefficients, we have 
B—6o2 = 5A,,  208-+80* = 5Ay, (19) 
B?— 16078 = dy, 4af? = x. (20) 
From (19) we get 40°+-2A,a—A, = 0 (21) 
and B = 6a7+5A,. (22 
Substituting for 8 in (20) and reducing by means of (21), we have 
x = 36A,07+16A$a+42A, A, | (23) 
y = 2A, 02—3A,a+5A3 J 


The three values of «, real or imaginary, given by (21) determine, 
when substituted in (23), three real or imaginary nodes. Pliicker’s 
equations confirm that these are the only double points. 

By (17), (22) we have for the pair of nodal parameters 


?—2ot+6a?+5A, = 0. (24) 
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These are real if and only if 

a?+A, < 0. 
But by squaring and rearranging in (21) we get 

(a?-++A,)(16a4-+4A3) = 443+ A3. 

Thus the tangents are real at a real node, if and only if 

4A3+A3 < 0, 
which is the inequality (16) already obtained as the condition for 
three real cusps. Actually the limiting case 443+ A} = 0 gives co- 
incident tangents at the node, as it also gives coincident cusps. As 
we shall see, in this case two cusps and a node have moved up to 
coincidence in a multiple point of the envelope. We further remark 
that the roots of (21) will be all real if 

8A3+27A? < 0, (26) 

an inequality necessitating (16) but not necessitated by it. Of these 
two inequalities it is (16) which is the important one, as we shall see 
shortly. We have proved, then, the following: 

A. (i) if 443+A? < 0, there are three real cusps and one or three 
real nodes each with real tangents; 

(ii) of 443+ Az > 0, there is one real cusp and one isolated point, 
but no other real singularity. 

As an additional check on the figure, it is helpful to know the 
nature of the tangents (other than the cuspidal tangents themselves) 
that can be drawn from the cusps to the curve. If the cusp is at 
¢ = 7 and we write 7 for ¢ in (14), (15), we can substitute for x, y, A, 
in the quintic, getting, since a cusp represents a triple root, the triple 
factor (t—7)*. By consideration of coefficients the residual factor is 


found to be #24 3rt-+672-+ 10A,, (27) 


which represents the two tangents that can be drawn from the cusp 
to some other part of the curve. They are real, if and only if 


377+8A, < 0. (28) 
Writing 7 for ¢ in (15) and squaring, we get 
7§+6A,74+9A37? = Aj. 


Writing 37?+-8A, = —80, we get 
6+ 6A, 62—15A, 0+ (843+ 2743) = 0. 
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Now (i) if 8A3+27A? < 0 and therefore incidentally A, < 0, then 
by Descartes’ rule of signs we see that (29) has at most one positive 
and two negative roots 6; 

(ii) if SA3+-27A? > 0, the rule of signs gives at most two positive 
roots and at most one negative root 6. 

This argument is of interest only in the case of three real cusps, 
so that the above numbers of positive and of negative roots are 
actually attained. Recalling that 8A}+-27A} < 0, as part of the 
inequality (26), is the condition for three real and distinct nodal 
points, we may state our result thus: 

B. (i) if the curve have three real and distinct nodes, only one cusp 
is convex to the curve; 

(ii) if the curve have only one real node, two cusps (and no more) are 
convex to the curve. 


Forms of envelope: the conic 0 = 0 


With the foregoing considerations A and B to guide us, we find 
that the envelope has three possible types of figure: Figs. 3, 4, 5, 


Fic. 3. (443+ A? > 0). 


together with the two intermediate limiting forms Figs. 3-4 and 4-5, 

where, in Fig. 3-4, two cusps, a node, and the region [5] of Fig. 4 

have all shrunk to the single multiple point P, and in Fig. 4-5 two 

nodes, a cusp, and one of the regions [3] have shrunk to the single 

point Q. In the several figures the regions [1], [3], [5] have been 
3695.5 Co 
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determined by inspection, and we may note that those cusps which 
are described in B as being convex to the curve are, of course, the 
cusps which belong to the regions [5]. 

We observe that, in Figs. 4, 5, and 4-5, the regions [5] and [1] are 
mutually accessible by passage through the nodes of the envelope, 


Fic. 4. (443+ 42 < 0 < 84342742). 


and we consequently need, for the complete discrimination of the 
roots, some further curve of separation passing through the three 
nodes. The nodes are given by (21), (23), i.e. by 

x = 36Az07+ 1645 0+42A, Ag, 

y = 2A,0®—3A,0+543, 
where 403+ 2A,a—A, = 0. 
Any a-eliminant of these three equations gives a curve through the 
nodes. I choose the particular conic 

© = 2?+16A,y?—76A,A,2— 
—(272A3—108A2)y+24A2(4043+2742) = 0. (30) 
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Now it is necessary not only that our curve of separation pass through 

the nodes but also that it cut the envelope in no real points other 
than the nodes. Elimination of a between (21), (24) gives 

28 224, t4+- 144,04 68424724, A,t+40434+2742 = 0 (31) 

as the sextic whose roots are the six parameters of the three nodes. 

If we substitute for x, y from (14) into (30), we obtain an equation 


Fic. 5. (843+2742 < 0). 


in t of degree 10 which includes (31) and contains the residual factor 
8(#4—2A,t+3A42) = 0. (32) 
The discriminant of this quartic is a positive numerical multiple of 
16A8—A$ 


and is therefore positive if 
4A34+ Ai < 0. 
Moreover, its Hessian vanishes. Hence, by the conditions of (11) 
above, none of its roots are real, and © = 0 is therefore a suitable 
curve of separation. 
Other conics through the nodes are 


4y*— A, x—36A} y+80A$+54A, Aj = 0, (33) 
xy—4A3x—18A, A,y+40A} A,+27A} = 0, (34) 
a2—72A, A,a—(128A3—108A2)y+640A3+43242 A? = 0. (35) 
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The corresponding residual factors are, 
for (33), 2(#+A,) = 0; for (34), —2+A, = 0; 
for (35), 8(t4—A,?t?—2A,t+2A3) = 0. 
In none of these three, however, is the condition 443+ A? < 0 suffi- 
cient to secure all roots imaginary. The equation (32) which I have 
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chosen is a combination of (33), (35); it is possible that other more 
suitable combinations of (33), (34), (35) might be found, but I have 
not thought it worth while to spend time over the pursuit. 

The relevant branch of the hyperbola © = 0 is shown in broken 
line in Figs. 4, 5, and 4-5. The region [5] is on the concave side of 
this hyperbola, that is, on the side © < 0. We can read from the 
figures the following necessary and sufficient criteria for real and 
imaginary roots of the quintic: 

[5] Five real 

A>0, 448442 <0, 
A= 0, 4A3+ A> = 0, 
[3] Three real, two imaginary 
A< 0; 
A= 0, 4A3+ Az > 0; 
4A3+ A? = 0, 0 A~0; 
4A34+ Az < 0, 0>0;) 


————-, -—__—__ 


[1] One real, four imaginary 
A> 0, 4A3+ Az > 0; | 
4434+42<0, O>0.) 


(38) 


These alternatives between them exhaust all the possibilities. It is 
remarkable that the inequalities (36), (37), (38) are exactly equivalent 
to the corresponding inequalities (9), (10), (11) for the quartic. We 
have only to replace the A, and 2—9A3 of the earlier inequalities by 
443+ A? and © to obtain the inequalities for the quintic. 


A numerical example 
To illustrate the foregoing discussion of the general quintic we may 
apply the method briefly to a particular example in which A,, A; 
have convenient numerical values. I take the very simple instance of 
&— 108+ 5yt+a = 0, (39) 
in which A, = —1, A, = 0. The envelope is defined by 
a= 4(/—58), y= —(t—6P). 
The cusps, given by ¢ = 0, +3, are at (0,0) and (+ 24v3,9). The 
nodes, given by t= +v5 and +v2¢ = +v5+1, are at (0,5) and 
(+82,4). Taking eight z-units equal to one y-unit, we get the 
envelope roughly as in Fig. 6. 
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The hyperbola 0 = 0 is 
x?— 16y?+ 272y— 960 = 0. (40) 
But let us consider more generally the conic 
a2— 16y2+ 272y—960-+32\(y—4)(y—5) = 0, (41) 
which also passes through the three nodes, since the line-pair 
(y—4)(y—5) = 0 passes through them. Substituting for x, y in terms 


Fia. 6. 


of ¢t and removing the factor (¢?—5)(t4—6t?+-4) = 0 corresponding to 
the three nodes, we have the residual factor 
4+ 2r7+-2A+-3 = 0. 
The roots of this quartic are certainly imaginary if A” < 2A+-3, i.e. if 
—1<A<3. They are also evidently imaginary, if A > 0. Thus 
(41) provides a more general curve of separation, if A > —1. 
If we apply Sturm’s method to (39) and seek the condition for all 
roots real, we get (excluding limiting cases) besides A > 0 the two 


inequalities y¥ <6, «2— 16y?-+ 176y2—480y < 0 


instead of the single inequality © < 0. The second of these Sturm’s 
inequalities introduces a bipartite cubic passing through the cusps 
as well as the nodes (it is not difficult to see that this is always the 
case with Sturm’s penultimate remainder). This cubic is shown in 
Fig. 6 by the dotted line: the other Sturmian boundary y = 6, which 
is the tangent to the infinite branch of the cubic at its ‘vertex’, is 


also shown in dotted line. The conic © = 0 is shown in broken line. 
4 











GROUPS WHOSE SYLOW SUB-GROUPS OF ORDER / 
p" CONTAIN NO MORE THAN p”*? OPERATORS 
WHOSE ORDERS ARE POWERS OF p 


By G. A. MILLER (Urbana, I.) 
[Received 4 August 1933] 


Ir a group G involves a Sylow sub-group H of order p™ which has 
at least two different cross-cuts of index p under it with its other 
conjugates, then H must also have cross-cuts of index p* with at least 
p* of its other conjugates. This results directly from the fact that 
one of the p other conjugates of H which have a common cross-cut of 
index p under H with H must also involve two different sub-groups 
of index p under it which are cross-cuts with its other conjugates. 
One of these is not in H but has a sub-group of index p? under it 
which is contained in H. The transforms of this Sylow sub-group by 
the operators of H are p* other conjugates of H which have cross-cuts 
of index p? under H with H. This includes the theorem that if all 
cross-cuts of H with its other conjugates are of index p under H, 
they are identical. 

The given p? conjugates of H which have cross-cuts of index p* 
under H with H involve p sub-groups of order p”-! which contain 
these cross-cuts and are contained in the p given sub-groups of order 
p” which are conjugate under H and have a cross-cut of index p under 
H with H. If the given cross-cuts of index p* under H were not 
identical, they would generate this cross-cut of index p under H. 
These p® conjugates of H under H could then have no common 
operators besides those which appear in H, since such a common 
operator would transform into itself this cross-cut of index p under 
H, and hence a Sylow sub-group of order p™ in which it does not 
appear. As this is impossible, the following theorem has been estab- 
lished: A group which contains a Sylow sub-group H of order p™ which 
has at least two different cross-cuts of index p under it with its other 
conjugates contains more than p™*® operators whose orders are powers 
of p whenever a set of its p* conjugates of H which are conjugate under 
H does not have a common cross-cut with H of index p* under H. 

A proof of this theorem may be based upon the fact that the 
p”*2—pm*+! operators which were proved to be distinct in the pre- 
ceding paragraph do not include any of the p+! operators in the 
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given set of p-+ 1 Sylow sub-groups including H which have a common 
cross-cut with H of index p under H, and this set involves exactly 
p”*! distinct operators. These p™+?—p™+! operators do also not 
include any operator except those of H in the other set of p+1 Sylow 
sub-groups including H which have a different common cross-cut 
with H of index p under H. The last statement results from the fact 
that the given set of p? Sylow sub-groups does not involve the cross- 
cut of the two given different sub-groups of index p under H, since 
this cross-cut is invariant under H, while the cross-cuts of H with 
the given p*? Sylow sub-groups are not invariant under H since these 
p* sub-groups were assumed to have different cross-cuts with H and 
also to be conjugate under H. 

Continuing the assumption of the preceding paragraphs, except 
that the p cross-cuts of H with the given p? conjugates of H are now 
supposed to be identical, we proceed to consider the smallest possible 
number of operators whose orders are powers of p which can appear 
in these p* conjugates of H and in the set of 2p+-1 Sylow sub-groups 
of order p™ composed of H and the two given sets of p sub-groups 
which are conjugate under H but involve different cross-cuts with 
H of index p under H. We determine this smallest number by 
noting the number of the operators of each one of the given (p+1)? 
Sylow sub-groups of order p™ which cannot appear in any other such 
sub-group. This number is obviously p”—(2p—1)p”-*, since no more 
than two distinct sub-groups of index p of one of these Sylow sub- 
groups can appear in others of these (y-+1)* Sylow sub-groups in- 
cluded in these two sub-groups. Hence there results the following 
theorem: If a group involves at least two Sylow sub-groups of order p™ 
which have different cross-cuts of index p with one of their conjugates, 
then the group involves at least p"** operators whose orders are powers of p. 

Suppose that G involves at least one Sylow sub-group H of order 
p™ which has a cross-cut of index p under it with one of its conjugates, 
but that all its cross-cuts of this index with its conjugates are identi- 
cal. Since H transforms these conjugates among themselves, the 
number of the Sylow sub-groups of G which involve this cross-cut is 
of the form 1-+-kp. As H transforms the remaining Sylow sub-groups 
of order p” contained in G in sets of at least p?, it results that the 
number of these Sylow sub-groups in G is (1+kp)(1+/p?). No two 
operators of such a set of 1+kp Sylow sub-groups which do not 
appear in the same Sylow sub-group could generate a group whose 
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order is of the form p%, since this group would transform the common 
cross-cut of the set into itself and would therefore have a sub-group 
of index p in common with this cross-cut. 

We can now find as follows a lower limit of the number of the 
operators whose orders are powers of p which appear in such a group: 
Arrange the distinct sets of 1+kp Sylow sub-groups of order p™ 
linearly. The first set clearly contains kp™+!+(1—k)p™ distinct 
operators whose orders are powers of p. Every Sylow sub-group of 
the second set has at most p”-® operators in common with those 
of the first set, since all such commor operators must appear in the 
same Sylow sub-group of the first set. Hence the first two sets con- 


tain at least 2hp™+14 2(1—k)p™—(1+kp)p™-? 


distinct operators whose orders are powers of p. Similarly, the first 
three sets may be seen to contain at least 
3kp™+1+ 3(1—k)p™—3(1+-kp)p™-* 
distinct operators whose orders are powers of p, while the first n sets 
contain at least 
nkp™+1 + n(1—k)p™— 4n(n—1)(1+kp)p™ 
such distinct operators. When n = p+1, this becomes half of 
2hp™+2+ (2—k)p™+1+4 (1—3k)p™—p™-, 
and hence it is at least equal to p”*? when k is any natural number. 

It has now been proved that when G contains a Sylow sub-group 
H of order p™ which has a cross-cut of index p under it with some 
of its conjugates but not with all of them, then the number of the 
operators of G whose orders are powers of p must be at least equal 
to p™*?, and that when H has cross-cuts of index p under it with all 
of its 1+kp conjugates, then these cross-cuts are identical and the 
number of the operators of G whose orders are powers of p is of 
the form kp™+1+-(1—k)p™. 

We shall now consider the case when no two Sylow sub-groups of 
order p™ contained in G have a cross-cut which is of index p under 
them. It will be proved that in this case the number of the operators 
of G whose orders are powers of p is also at least as large as p™**. 

Let H be a Sylow sub-group of order p” contained in a given group 
G, and let p* (a < m—1) be the largest cross-cut of H with its other 
conjugates, while p*+? is the order of the normalizer of this cross-cut 
under H. The totality of the operators of G which transform this 
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cross-cut into itself constitutes a sub-group of G involving more than 
one Sylow sub-group of order p**?. The p™-*-8 conjugates of this 
sub-group of G under H have the property that no two of their Sylow 
sub-groups of order p*t® contain any common operators besides 
operators of H, since no two Sylow sub-groups of G have more than 
p* common operators. Each of these sets of conjugates contains p? 
Sylow sub-groups of order p**? besides the one found in H. Hence 
the total number of the distinct operators contained in all of these 
sets of conjugates but not found in H is p”(p8—1). When B > 1 the 
number of the operators of @ whose orders are powers of p must 
therefore be at least equal to p”*+*, and hence we may assuime in 
what follows that B = 1. 

When « = m—2, then H transforms one of its conjugates into p? 
others which have p distinct cross-cuts of order p* with H, since 
8 = 1. The number of the conjugates of H in G must exceed p?+-1, 
and hence it must be at least 2p?-+- 1, for, if it were p?+-1, a conjugate 
of H which differs from H would involve a sub-group of order p* 
which would not appear in H but would appear in p+1 of its con- 
jugates. This is impossible, since no two of these conjugates could 
appear in a set of p Sylow sub-groups of order p” which has a common 
cross-cut of order p* with H but does not include H. If the given 
set of p?+1 Sylow sub-groups of order p” would involve no common 
operators except those of H, the set would involve p™+? distinct 
operators, and hence it will be assumed in what follows that other 
common operators appear therein. The given p cross-cuts of order 
p* which are conjugate under H have themselves a cross-cut of order 
p”-%, which is invariant under H because it is invariant under two of 
these cross-cuts. These cross-cuts are invariant under each other. 

The said other common operators transform this cross-cut of order 
p”- into itself, since this is also common to the two given Sylow 
sub-groups. When H has a cross-cut of order p* with a Sylow sub- 
group of order p” which is not contained in the given set of 1+ p? 
such sub-groups, then it results that G involves more than p”+*? 
operators whose orders are powers of p, since 2p? Sylow sub-groups 
appear in sets of p such that each set has a cross-cut of order p* with 
H. Hence one of these Sylow sub-groups cannot have a cross-cut of 
order p* not contained in H with more than 2p—1 of these sub-groups 
besides itself. The common cross-cuts of order p”-* of the two sets 
of p* Sylow sub-groups of order p” may be assumed to be the same 








GROUPS OF A CERTAIN TYPE 27 
since otherwise the number of operators whose orders are powers of 
p in these two sets would obviously exceed p”+?. Hence no more 
than (2p—1)(p”™-*—p™-*) of the operators of a Sylow sub-group 
which are not found in H could appear in others of these 2p? Sylow 
sub-groups. The number of the operators of G whose orders are 
powers of p must therefore be at least 

2pet)— 2p ip— i(p—1)—9* = Bg — ig + oe a, 
and hence at least p™*+?. 

- When H has no cross-cut of order p* with a Sylow sub-group of 
order p™ which is not contained among the given set of 1+p? such 
sub-groups, it has the given invariant sub-group of order p”-* as 
a cross-cut with p* other conjugates, since one of the conjugates of 
H which differs from H and appears in the given set of p?+-1 con- 
jugates of H has a cross-cut of index p*? which involves this invariant 
sub-group of order p”-* with a sub-group of order p™ which is not 
found among these p?+-1 such sub-groups. Moreover, the cross-cuts of 
H which appear in these p?+-1 sub-groups generate a group of order 
p”™— under which this sub-group of order p™ has p? conjugates com- 
posed of those which have cross-cuts of order p”-* with the groups 
in the given set of p?+1. As these sub-groups have only the given 
cross-cut of order p”-* in common with the rest of the said p* con- 
jugates, it results as above that @ involves at least p™**® operators 
whose orders are powers of p. 

[t remains to consider the case when a < m—2. In this case we 
may consider p?+1 conjugates of H composed of H and p* other 
Sylow sub-groups of order p” which have cross-cuts of order p* with 
H and are conjugate under a sub-group of order p*** contained in H. 
These appear in p? sets of p each such that each set has the same 
cross-cut of order p* with H. If one of these does not have a cross-cut 
of order p* which does not appear in H with another of the set of p* 
Sylow sub-groups, then the set clearly contains more than p”*? distinct 
operators whose orders are powers of p. If it has such cross-cuts, then 
it cannot have more than p?—1 of them, since it cannot have such a 
cross-cut with more than one in a given set of p having the same cross- 
cut of order p* with H. It can therefore not have more than p™-* 
of its operators which do not appear in H in common with its other 
conjugates, since its cross-cut with one of these conjugates cannot be 
of a larger order than p”-!. This completes a proof of the following 
theorem: If the order of the largest cross-cut of a Sylow sub-group of 
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order p” with one of its conjugates is less than p™-|, then the group con- 
tains more than p™** operators whose orders are powers of p. 

Suppose that G involves exactly (1+-p)? Sylow sub-groups of order 
p™ which do not have a common cross-cut of index p under each of 
them. Hence H involves at least one sub-group of index p? under H 
which is the cross-cut of H with at least p of its other conjugates, 
and it also involves two sub-groups of index p under it which are 
separately the cross-cuts of H with p of its conjugates. These two 
sub-groups must be distinct, since p” is not divisible by 1+2p. The 
cross-cut of these two distinct sub-groups is of order p”-® and in- 
variant under H. Each of the 2p Sylow sub-groups with which H has 
cross-cuts of index p under H has p cross-cuts of order p”—! with p of 
the p? Sylow sub-groups of order »” which are conjugate under H. 
Hence the given cross-cut of order p”~* is the common cross-cut of H 
with the p? Sylow sub-groups of order p” which are conjugate under H. 

It is now easy to determine the number of the operators whose 
orders are powers of p which appear in such a G. The given sub- 
group of order p”~? is invariant under G since it is the maximal 
cross-cut of all the Sylow sub-groups of order p” contained in G. 
Hence each of the (p+1)? Sylow sub-groups of order p” contains 


p™-*(p—1)? operators which appear in no other Sylow sub-group of 
this order. It also contains 2p”-!—2p”-* operators which appear in 
exactly p-+-1 such Sylow sub-groups, while its remaining p”-? opera- 
tors appear in each of these Sylow sub-groups. This proves the 
following theorem: Jf a group contains exactly (p+1)? Sylow sub- 
groups of order p™, and if not all of these sub-groups have a common 
largest cross-cut, then the group contains exactly p™** operators whose 


orders are powers of p. 

From this theorem it results directly that the number of the opera- 
tors whose orders are powers of p which appear in a group does not 
necessarily increase with the increase in the number of its Sylow 
sub-groups of order p”, since a group which has exactly p?+1 Sylow 
sub-groups of order p” which have a common cross-cut of order p”~? 
with one of these sub-groups obviously contains also exactly p™** 
operators whose orders are powers of p. For instance, the icosahedral 
group contains exactly 16 operators whose order are powers of 2, 
while the direct product of two non-cyclic groups of order 6 has the 
same property. The former contains 5 Sylow sub-groups of order 4, 
while the latter contains 9 such sub-groups. 
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It was noted above that when a Sylow sub-group of order p™ has 
cross-cuts of index p with all of its conjugates, then these cross-cuts 
are identical and this common cross-cut is an invariant sub-group of 
the group G. Hence it results that if every Sylow sub-group of G has 
a cross-cut of prime index with all its other conjugates, then G con- 
tains an invariant sub-group, which may reduce to the identity, 
corresponding to each of its Sylow sub-groups. Every two such 
invariant sub-groups have only the identity in common, and hence 
G involves as an invariant sub-group the direct of all these invariant 
sub-groups. The order of the corresponding quotient group is then 
not divisible by the square of a prime number, and hence G must be 
solvable. From this there results the following theorem: If every 
Sylow sub-group of a group has a cross-cut of prime index with each 
of its conjugates, then the group is solvable and involves an invariant 
sub-group which is the direct product of sub-groups of prime index under 
each of its Sylow sub-groups. It should be noted that this theorem 
includes the theorem due to G. Frobenius that a group is solvable if 
its order is not divisible by the square of a prime number. 
A corollary of this theorem may be stated as follows: If a group 
contains less than p™** operators whose orders are powers of p for every 
prime number p which divides its order, then the group is solvable, p™ 


being the highest power of p which divides its order. In some of the 
cases considered above, the number of operators whose orders are 
powers of p was obtained by noting such operators which appear in 
only one of the Sylow sub-groups considered. Hence it should be 


noted here that a Sylow sub-group does not necessarily involve any 
operators which do not appear in any of its other conjugates. For 
instance, in the symmetric group of degree p? (yp > 2) all the substitu- 
tions of any given Sylow sub-group of order p”*+! appear in the other 
Sylow sub-groups of this order since such a Sylow sub-group may be 
transformed by a permutation of order 2 which transforms a cycle 
of order p into its inverse and involves only letters of this cycle. Such 
a permutation transforms into itself the direct product of p regular 
sub-groups of order p but does not transform the Sylow sub-group 
into itself.* Since each of the operators of this direct product appears 
also in other Sylow sub-groups, it results that there are groups in 
which each of the operators of a given Sylow sub-group appears in 
some other Sylow sub-group but not in the same other such sub-group. 
* G. A. Miller, American J. of Math. 23 (1901), 176. 











‘A ONE-DIMENSIONAL UNIVERSE OF DISCRETE | 
PARTICLES 


By E. A. MILNE (Ozford) 


[Received 13 October 1933] 


In a recent paper* I have shown that, if a continuous one-dimensional 

distribution of moving particles is to have the same description as 

viewed from any particle of the system, then the number of particles 

with velocities between u and u-+-du, reckoned with respect to any 
€ ~ 7 Ss 7s / m Ta 

particle of the system, is i ie 


c(1—u?/c?)” 
For the functional equation 


f(w’) dw’ = f(u) du, 
‘= u—t - du’ = sei te ale 
1—uV /c? (l—uV /c?)? 
and V is arbitrary, has no other solution independent of V; (2) is the 
condition of identity of description of the velocity-distribution as 
described by two particle-observers whose relative velocity is V, and 
(3) is Einstein’s law of composition of velocities. 

The object of this note is to obtain the corresponding distribution 
for a set of discrete particles in uniform motion. The interest of the 
result is that it gives insight into the structure of that idealized 
representation of the universe in which the extra-galactic nebulae are 
replaced by discrete particles in some three-dimensional distribution. 
The problem of the continuous three-dimensional distribution of 
moving particles such that the distribution is described in the same 
way by any particle-observer of the system was also solved in the 


where U 


paper cited, the solution being shown to be 
A dudvdw 
ra RT 
u, v, w being the Cartesian components of velocity. I have not suc- 
ceeded in solving the problem of the discrete distribution in three 
dimensions, and I do not know whether it possesses a solution. If 
it does possess a solution, we should thus obtain a very valuable 
world-model. If it does not possess a solution, highly interesting 


* Zeits. fiir Astrophys. 6 (1933), 35. 








A ONE-DIMENSIONAL UNIVERSE 31 
consequences would follow, for the description of the velocity- and 
space-distributions would then differ from particle to particle, and 
the residual differences would be expected to be associated with resi- 
dual gravitational fields,* and the tidal and other influences of these 
might well be responsible for the occurrence of rotation, the existence 
of a plane of symmetry, and the formation of spiral arms (as sug- 
gested by Jeans) in the extra-galactic systems. 

In default of a solution of the three-dimensional problem for dis- 
crete particles, the solution of the one-dimensional problem is not 
without interest, and it is of interest also to verify that the discrete 
distribution passes continuously into (1) as the particles crowd 
together. 

Suppose, then, that as recorded by a given observer, attached to 
a given particle, the velocity-distribution is specified by the sequence 

— eee ee Se Se ee 
where uw) = 0 is the observer’s velocity and the w’s are in ascending 
order of magnitude. 

Shift the observer to the particle moving with velocity u,. Then 
the particle previously observed as moving with velocity wu, ,, is now 
observed to be moving with velocity wj,,,, where by (3) 

Unsi— “4 

1—; Uy 43/C?” (*) 
But if the velocity-distribution is to be described in the same way 
by the new observer, we must have 

Un+1 = Uy» (5) 
for the nth particle in the new view must have the same velocity 
relative to the new observer as the nth particle in the old view had 
to the old observer. Moreover, if condition (5) is satisfied for all n, 
then the condition of identical description will be satisfied for any 
shift of the observer’s head-quarters. For if he shifts from uw, to u,, 
this shift can be accomplished in the stages uw. to uw ,, Uz tO Ug,..., sy 
to u,, and each of these shifts leaves the description of the velocity- 
distribution unaltered, for each shift is a reproduction of the first 
shift. Hence (5) is a sufficient condition. Inserting it in (4), we have 


Uns a 


Unsi—4 (6) 


uu, = . 
2 
1—Uy Uy, 44/€ 


n 


* The residual field is, of course, zero at each particle when the problem 
has a solution. 
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Put «,,/¢ = a,. 


= * | le 
Then (6) becomes “a, 
i— My Xney 


which gives, on solving for «,,,,, 


| 

: Xn Ty 

Xn +1 ~ l 3 . 
s Xy Xn 


Hence ss (1+o,)(1-+a,) 


I Ty aX, 


(l—a,)(1—«,) 


> 


1+ Xy X, 
ee l+oa l+-a, l+a 
or, dividing, are eee B nn, 
Saat eee — Xn l—a, 


By repeated application of this relation we find 
Ions (7 yn" 
~ > 
I—a45 l—a, 
1+-u,/c)"—(1—u,/c)” 
whence = of Be ” : a (10) 
(1+ 4,/c)"+(1—wu,/c) 
This determines the invariant one-dimensional discrete velocity- 
distribution w, in terms of a single arbitrary constant u,. The central 
velocities are given by 
U_» = —2u,/(1+-u?/c?), U_, = —th, 
Mu, = 0, Uy = 2u,/(1+-u§/c*), Uz = U,(3+u3/c?)/(1+3u2/c?). 
As n>, u,—>c and u_,—>-—c. For u,/c small and n finite, 


Uy, ™~ NU. 


a 


To compare this distribution with the corresponding continuous 
distribution we proceed thus. The index n of the particle moving 
is given by 


n= (log V(ee; — ‘), 
—U,, | 1—u,/c 


whence, when the particles crowd together in velocity and successive 


with velocity u 


a 


u,,s are close together, 


| lp\ —1 
dn a 1 (lee ee) (11) 


S ; 2 /p2 . 
ou, cil—u;/c 1—u,/ 


The number of particles in the range of velocity (u,u+dwu) is then 
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(5n/5u) du, where we can put u for uw, in (11). Thus it is 
. a _. a (12) 
log(1-+-,/c)—log(1—1u,/c) c(1—u?/c?) 
This agrees in form with the continuous solution (1), and shows, 
moreover, that l 
sie log(1-+-u,/c)—log(1—u,/c) 
It follows that the discrete distribution mimics a continuous distribu- 
tion (a) for any u,, when » is large; (b) for any uw, when w, is small. 
The foregoing algebra suggests that if a three-dimensional invariant 
discrete distribution exists, its description will involve three arbitrary 


vectors uy), u'?), u‘>), 
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’ ON THE ASYMPTOTIC PERIODS OF MERO- 
MORPHIC FUNCTIONS 
By J. M. WHITTAKER (Liverpool) 
[Received 8 November 1933] 


1. An asymptotic period of an integral or meromorphic function f(z) 
was defined in a recent paper* to be a number f (+ 0) such that 
Ag(z) = f(z+B8)—f(2) (1.1) 
is of lower order than f(z). A familiar example is the Weierstrass 
function ¢(z), for which 
C(z+2u,)—L(z) = 2m, — E(e-+2ug)—L(z) = 2mq, (1.2) 


so that the numbers 2kw,+2/w,, where k, | are any integers, are 


asymptotic periods. 


The following result was established: 

THEOREM A. The asymptotic periods of an integral function lie on 
a straight line through the origin and form a set of measure zero. 

A partial extension of this result to meromorphic functions is 
proved below. 

THEOREM B. The asymptotic periods of a meromorphic function 
either 

(i) form a denumerable set; 


(ii) lie on a straight line through the origin. 

In case (i) it is shown by means of examples that the set may be 
everywhere dense in the plane or everywhere dense on each of a set 
of equidistant parallel lines. 

In case (ii) the set may be non-denumerable.+ Theorem A suggests 
that it must have measure zero, but this I am unable to prove. 

The proof of Theorem A consisted of six stages, the first concerned 
with proving that the set must lie on a straight line through the 
origin, the others with proving that it must have measure zero. When 
it is attempted to apply the same method to meromorphic functions 

* Whittaker (5). If f(z) is of infinite order, B is an asymptotic period when 
A®(z) is of finite order. 

+ An example of an integral function for which this is the case is given in 
Whittaker (6). To obtain an example of a meromorphic function it is only 


necessary to add z7}. 
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serious difficulties of principle are encountered and, except in the 
first stage, I have not been able to overcome them. 

Theorem B is the conclusion of this first stage. It is not nearly so 
‘deep’ as Theorem A, but it does at least show that the asymptotic 
periods of a meromorphic function always form a set of plane measure 
Zero. 

2. The proof depends on a ‘geometrical’ argument based on the 
following lemma. 


LemMMA 1. Let 4(z) be an integral function of finite order p. Then 
corresponding to each number o > p there is a number d, such that 


ff log+ |d(z)|-1dS < 3 r+ D°8(A) (d> d,), (2.1) 


where A is any domain in the complex plane, S(A) its area, r,, r9,... the 
moduli of the zeros of 4(z), the summation is taken over the zeros inside 
or on the boundary of A (multiple zeros being counted multiply), and 
d = min|z|, D = max|z| (z in A). 
The function ¢4(z) is of the form 
P(zjexp(k,2+k,2z?+...+, 2°), 
where P(z) is a canonical product and s < p; so that, writing |z| = r, 
log+ |¢(z)|-? < log+| P(z)|-*+- |k,|r+...+ |k, |r*. 

Now 

[{ (iky\r+...+ [kelr8) dS < 4D*S(A) 


and* 
Ny Ss 
log|P(z)| > —KI+log | | dae 
1 
where k > 1. Also 


- m(a) rP + 


2B dz < rPt« (r 
xPtlytr 


: 


p being the genus of P(z). But (2.3), (2.4) give 
N N 
logt|P(z)|_ < KI+ } log+r,+ > log* 
1 1 
- 
< $4D°+ 3} log+|z—a, |“ 
1 


* Valiron (2), 53. 
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and so 
( log+|P(z)|1 dS < 4D°S(A)+ >| [ log+|z—a, |-?dS. (2.5) 
‘A A 
Now surround each zero a, with a circle C, of radius rz" (h > 4c), 
and let A,, be the domain A with C, excluded. Then 


( [ logt|z—a,,|-1 dS < hlogr, S(A) 
An 
and, if r, > 1, 

ryt on 


J2dS = | | logus udude 


Cr 0 0 


( [ logt |z—a 


Nr 
= n(hlogr, —43)r,™ 


“re (d > d>). 


On combining this inequality with (2.2) and (2.5) we get the result 
stated. 

We need also a lemma which assigns a lower limit to the ‘surface 
density’ of log+|¢(z)|. This was mentioned incidentally in some work 
in which more precise results of the same kind played an important 
part.* 


LemMa 2. Let 4(z) be an integral function of order p and let 


] 


S(a,r) = n(r?—a®) 


+|h(ue") \u dud. 


Then ee) 2 
roe logr 

3. Theorem B can now be stated in a more precise form. 

THEOREM. The asymptotic periods of a meromorphic function f(z) 
form a set B with one of the following properties: 

(ax) Bis null; 

(8) B consists of a set of points kw, k = +1, +2....; 

(y) B lies on a straight line passing through the origin and is every- 
where dense on that line; 

(5) B consists of a set of points kw+-lw’ (k,l = 0, +1, +2,...) where 
w’/w is not real and the case k = 0, 1 = 0 is excluded; 


* Whittaker (7), Lemma 6. Other results of the same type are given in 


(3), (4). 
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(c) B lies on a set of lines y = mx+kc, where m, c are fixed and 
k = 0, +1, +2,..., and is everywhere dense on each of them; 

(C) B is everywhere dense in the plane. 

Let p (possibly 0) denote the order of f(z) and « the exponent of con- 
vergence of its poles,* so that x < p. Then 

(i) of « = p, Bis a denumerable set; 

(ii) of x < p, the only possible cases are (x), (8), (y), but in case (y) 
B is not necessarily denumerable. 

Any set of points in the complex plane with the property that 
‘if wu, v are any members of it, then —w, u+-v are also members’ must 
either have one of the properties (a),..., (¢) enumerated above [with 
addition of the point 0 in cases (8), (8)], or else consist of the point 0. 
For if 2 is such a set containing a point wu other than 0, consider the 
projections of the points of F on the line joining 0, uw and on the line 
through 0 perpendicular to this. It is easy to see that the projections 
have the additive property and so, making use of a result of a pre- 
vious paper,t the projections on each line must have one or other of 
the properties (a), (8), (vy). The different possibilities correspond to 
cases (a) (£), e.g. (8).on one line and (y) on the other gives (e), (y) 
on each line gives (f). It is evident that if to B there is added the 
point 0 a set with the required additive property is obtained. The 
first part of the theorem has therefore been established. 

Passing to the second part, take first case (i). The order of Ag(z) 
is not inferior to the exponent of convergence of its poles. Now these 
consist of the poles 6,, b,,... of f(z) together with the poles b,—8, 
b,—B8,... of f(z+f), and so have exponent of convergence p, unless 
some of the poles cancel out. Thus 8 cannot be an asymptotic period 
unless it is of the form 6, —b,. It follows that B must be denumerable. 
[t should be noted that not every number of the form b,—6, need 
be an asymptotic period. 

Now take case (ii). f(z) is of the form 

fie) = ©), 
g() 
where g(z) is a canonical product of order x < p and h(z) is an integral 
function of order p. If f is an asymptotic period, 
h(z-+B)g(2)—g(2+B)h(2) 
g(z)g(2+B) 


* i.e. the lower bound of numbers k such that > |b,|* converges, b,, by,... 
being the poles of f(z). + Whittaker (5), Lemma 7. 
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is of lower order than f(z), and so 
F(z) = h(z+B)g(z)—g(e+B)h(z) 
is of order A < p. 
Suppose, if possible, that y is an asymptotic period which does not 
lie on the line joining 0,8. Then 
G(z) = h(z+y)9(z)—g(z+y)A(z) 
is of order p < p. 
Let o be chosen so that 
K,A,u<a<p, (3.1) 
and let d, be the associated number when ¢(z) = g(z) in Lemma 1. 
Moreover, let D, > d, be chosen so that 
lg(z)|, |F(2)|, |G@(z)| <<" — (r > Dy). 
The identity , 
1 ‘—1 
f(@o+kB+ly) = f(%)+ 2 Aor sy)+ 2 Apleotly + tf) 


shows that, if D > D, and z is a point in the annulus [ 


(Dy < |z| < D), 
then 


f(@) = f%)+ 2 G(z.)iG(2.)9(%s+y)} A+ +2 P(e )g(%)g(%+B)}" 


where z, is in a ‘period parallelogram’ or cell tommeet with periods 
B, y) intersected by the circle |z| = D, and 2,, Zo,..., 2], 23... are points 


' ' 


‘conjugate’ to z in different cells in the annulus. The number of 
these points does not exceed Kg, D, where Kg, depends only on , y. 


Hence AN , p? Lofo”\ | —2 
Ifle)| < |fleo)| + Kp. De?” max|g(z") |, 
where z/, z;,... are all the points ‘conjugate’ to z in the annulus. 
On making use of the inequalities* 


q 
log* (ay Og +++ Og) < bi log*a,, 
1 


log+(S a) < > log*a,-+log q; 
this gives 
log* |h(z)| < log*|g(z)|+-log*| f(z)| 
< logt|g(z)|+log*|h(zo) | +-log* |g(zo) |“? +log Kg + 
+log D+ D°+-2 log+ max |g(z")|-1+-log 2. (3.2) 


* Nevanlinna (1), 14. 
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Let Z denote the point (or one of the points) z? for which |g(z7)|- 
is a maximum. As z moves about in a cell, Z will jump from one 
cell to another, and as z moves over a whole cell, Z will move over 
a number of disconnected domains A,, Ag,..., A, in different cells, 
such that 
S(A,)+S(A,)+...+S8(A,) = area of a cell = S. 
A single A, may consist of several disconnected parts in the same cell. 
Now log*[h(@o)| < Di < 4D" (D> Dy), 
so that (3.1), (3.2) give 
log* |h(z)| < 3log*|g(Z)|*+2D" (D>D)), 


and, if C is any cell in the annulus I, 


ff logt|h(z)| dS < 3 | i] log*|g(Z)|—1 dS +2D°8. 
C s=1 6 


On applying Lemma I, this gives (writing |b,,| = r,,) 
na k 
| | log*|h(z)| dS <3 > > rn°+D°S(A,)} +2D°S 
“C e=1* 4, 
< 123 r°+ KD? < KD», (3.3) 
n=1 


the factor 12 being arrived at by observing that a pole b,, may be 
on the boundaries of four different domains A,, if it happens to be 
the corner of a cell. From (3.3), 


[ { log*|h(z)| dS < KDen(D®—D3) (D> D,), 
J 


and so 
1 
lim (log D)-*1 log*+{h(z)| dS < ; 
r 
and this contradicts Lemma 2. Hence the assumption that y is not 
collinear with 0 and is false, and so either B is null or B lies on 
a straight line through the origin. 


4. Examples. An integral function cannot have an asymptotic 
period (and a fortiori cannot have a true period) unless its order is 
at least one; and a meromorphic function cannot be singly periodic 
unless its order is at least one nor doubly periodic unless its order is 
at least two. The examples which follow show that there is no such 
limitation on the asymptotic periods of meromorphic functions. All 
the cases (a),..., (¢) oecur among functions of any positive order. The 
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examples given are all of functions of order less than one, but it is 
easy to construct similar examples of higher order. 

x) If f(z) is any integral function of order less than one, then 
f(z)+2" is a meromorphic function with no asymptotic periods. 


(8) Write on 
$,(2) = | | (1555) (4.1) 


k 
where c = 24? (0< p< 1). Then 


i 2) 
$(z) =] ¢,(2) 
n=1 
is an integral function of order p and 


x(z) $(2) 


is a meromorphic function of order p. Now 


—d, ¢,(z+1) 
2)= —log 
Lan $,(2) 


4) 
( ] 
= > — ? 
cr Qn c”®—z 
n=1 r 


and this function is of order zero. Thus all positive and negative 
integers are asymptotic periods of x(z), and it is easy to see that they 
are the only asymptotic periods. 

(y) The remaining examples are constructed on the same principle, 
though the details are more complicated. If x(z) is defined by 


(4.2), (4.3), where now 
e; = ad - (4.4) 
»2n k 
On Ton 


it will be found that x(z) is of order p and that, if p is any integer, 


1 ’ 
x(e+ 35} xe) is of order }p. 


Thus the asymptotic periods are everywhere dense on the real axis. 

There are no complex asymptotic periods, for, since x = p, an asymp- 

totic period must be the difference of two poles, and these are all real. 
(5) If we replace (4.1) by 


9=T1( ayer) in: 
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then x(z+1)—x(z) is of order }p, 


x(z+7)—x(z) is of order }p, 
Thus the asymptotic periods are of the form k+-li, where k, 1 are any 
integers. 
(c) If we replace (4.1) by 


4” 
= 2 

#.(e) = | | (1-—4—} 

rea 4 ; 

k,t=1 c rte 

then, if p is any integer, 
] , 

x(: +35} —x@) is of order 3p, 
x(z+7)—x(z) is of order $p. 


The asymptotic periods are everywhere dense on each of the lines 


y = k (k = 0, +1, +2....). 
(¢) If we replace (4.1) by 


4” a 
o— 1 gex 
k,l=1 oe + oat on 


then, if p is any integer, 


x(= 4 )-x (z) is of order #p, 


x (2+ 35) —x@) is of order #p. (4.8) 


The asymptotic periods are everywhere dense in the plane. 
As a specimen of the proofs of these statements let us take (4.8). 
If p is any positive integer, 


p 


| 1 
+ 3) —x@) = > 5: [low da(2-+ 55) —lon but) + 


n=1 


where 


n=p+1 l=1 k=0 


“= TT] 


nm=p+1 l=1 k=4"—2"-?+1 
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and, as y,(z), w,(z) are integral functions of order 3p, this gives the 
result stated. 

5. There remains the problem mentioned in the introduction. Jn 
case (ii) is the set B necessarily of measure zero? It is easy to deduce 
from Theorem A that this is so if f(z) has only a finite number of 
poles. We need a further lemma. 

Lemma 3. Let h(z) be an integral function and g(z) a polynomial. 
Then there is an integral function Q(z) and a polynomial R(z) such that 

h(z) = Qz)g(z)+ R(z). (5.1) 

If g(z) = A(z—a)(z—B)...(z—v), 
then h(z) = h(«)+(z—a)h,(z), 
where h,(z) is an integrai function; 

hy(z) = hy(B)+(z—B)ho(2), 
and so on. A combination of these equations gives (5.1). To apply 
it, write a Q(z) 4.2, 
9(2) 
Then the asymptotic periods of f(z) are the same as those of the 
integral function Q(z) and so form a set of measure zero. 

The argument breaks down if the number of poles is infinite. For 
the essential point in it is that R(z) is of lower order than Q(z); and 
it is not necessarily true that, if h(z) is an integral function of order p 
and g(z) an integral function of order less than p, integral functions 
Q(z), R(z), the latter of order less than p, can be found satisfying (5.1). 


For example, let h(z) = e* 


g(z) = sin 72, 
and let Q(z), R(z) be any integral functions satisfying (5.1). Then 
R(n) = e”, 
so that the order of R(z) cannot be less than two. 
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MEAN-VALUE THEOREMS FOR THE 7 
RIEMANN ZETA-FUNCTION 
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[Received 3 January 1934] 


1. We consider in the present paper the existence of mean values of 
the form T 
lim 55 7 | me (o-+it)} dt = Melo) (o>}), (1) 
oe 
where {(s) = {(o+it) is the Riemann zeta-function. If o > 1, the 
existence of (1) for any function ®(z) of the complex variable z, con- 
tinuous except perhaps for z = 0, is very easy to prove. In what 
follows we obtain the existence of (1) also for } < o < 1 for certain 
(z).* 
The mean values which have been most discussed are mean values 
of the form , 
lim as [ \f(o-+-it) |?* dt = O(c). 


Ta = 
T 


Results concerning these mean values for integral values of k have 
long been known. Quite recently a new and general result concerning 
also non-integral values of k was given by Ingham.t 

In this paper we do not establish the existence of (2) for any new 
values of k and o, but we prove as a special case of a more general 


theorem that, if r 


] ee 
lim Vi | \f(a-+it) |?" dt < a (3) 
T-2 + 
A 
for a given o > } and an m > 0, then (2) exists for the same value 


of o and any positive k < m. 


* If the value of o is 1, so that the pole lies on the line of integration, 
it is always understood that an interval about 0, for instance the interval 

l t < 1, is left out by the integration. Once this remark is made we shall 
not repeat it. If (z) satisfies the condition ®(z) = ®(z), we can replace the 
mean over —T' < t < T by the mean over | < t < 7, but for general ®(z) 
these means are different. It could, however, be proved that in the cases in 
which we prove the existence of (1) we should have found the same value if 
the mean had been taken over 1 < ¢t < T. 

+ Ingham (3) (see the list at the end of this paper). The result generalizes 
some theorems of Titchmarsh ; for details and references see Ingham’s paper. 
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Our theorem is obtained as an almost immediate application of 
a previous theorem of the authors, given below as Theorem A, which 
states the existence of a distribution of the values of f(s) on any line 
s = o-+it, where o > 3; the mean values, of which we prove the 


existence, are expressed by this distribution. Ingham considers uni- 


formity in o in the existence of (2); we could do this, but should then 


first have to extend our previous result in that direction. 


2. We base our exposition on the following theorem.* 


THEOREM A. If o > } is fixed, then there exists in the complex plane 
of z = u+iv a continuous function F(z) = F(z;c) with the following 
property: if R is an arbitrary rectangle (u, <u < Us, Vy << VU < Vg) in 
the z-plane, with sides parallel to the coordinate axes, and if L(T') denotes 
the total length of those intervals in —T <t < T in which C(c+it) 
belongs to R, then 

MP) _ ifr 
= | | F(z) dudv. (4) 
> “- 

This theorem suggests the existence of the mean value (1) when- 
ever the integral of the function |®(z)| F(z) over the z-plane is finite, 
in which case we may expect to find 

Mso(c) = [ (2) F(z;0) dudv; (5) 


oa 


but we cannot conclude this from Theorem A except when {(o+7#) 
is bounded, that is, when o > 1. We shall, however, prove that the 
formula (5) holds in the cases in which we prove the existence 
of (1). 

Before we begin the deductions from Theorem A we observe that 
we must have 


[ F(z) dudv = 1. (6) 


* Bohr-Jessen (1), (2). Theorem A was not actually proved in these papers, 
which dealt not with the function €(s) itself but with log {(s), but it was 
observed that it could be deduced from the corresponding theorem for log &(s), 
which was proved in all details. It should be noted that a much weaker 

: ‘ ‘ ings ee ae 
theorem than Theorem A stating only the existence of the limit lim oT 
T+» 2 
for an everywhere dense set of rectangles R would be sufficient for the following 
deductions. This theorem is much easier to prove. 
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If not, we should have 
[{ Pe) dudv <a < l 
R 


for all rectangles R. Take R such that it contains the circle |z| < M; 
then it follows from (4) that we have, for any m > 0, 
- 
] 
| gf\|2 Pe 2m 

—-T 
for all sufficiently large values of 7’, which contradicts the fact that 
(3) is true for some value of m, e.g. for m = 1. 

3. We shall prove two theorems of which the first is an immediate 
deduction from Theorem A and the second (which is the main result) 
is a corollary of the first. 

THrorEM 1. Let o > } be fixed, and let D(z) be a continuous func- 
tion of z with the property that to any « > O corresponds an n = y(€) > 0 
such that, if 4(t) is a step-function which takes only the two values 0 and 
| and for which 


r 
lim Z [ p(t) dt < », (7) 
r+027' | 

—T 


] 


Tw 2 


- 
then we have lim Vi |D{l(o+-it)}|p(t) dt < «. (8) 
. 


Then the mean value (1) exists and its value is given by (5). 
Proof. Since (6) is true, we can choose, corresponding to a given 
0, a rectangle R such that 
[ [ F(z) dudv > 1—n, 
R 
where 7 = »(e) is the number given by the theorem. If we now take 
A(t) = 0 when {(o+it) belongs to R and ¢(t) = 1 elsewhere, the con- 
dition (7) is clearly satisfied. Hence, if we can prove that 


7 @ 


—T R 
for any continuous function ’(z), the proof will be finished. For 
then, taking first Y'(z) = |®(z)|, we infer the existence of the integral 


rs 
[ Vo+iyi—gey a= [[ VEFe) dud (9) 


[J ee) dude, 


—O 
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and then, taking V’(z) = ®(z), we conclude that (1) exists and is 
given by the formula (5). 

To prove (9) we consider, for some 5 > 0, a division of R into a 
finite number of rectangles R; of diameter < 5;* let A = A(d) denote 
the upper bound of |(z)—‘(z’)| for points z and 2’ satisfying 
lz—z’| < 6 and belonging to R; then, if z; denotes an arbitrary point 
in R,, the number 


> Fei) [ [ F "(z) dudv (10) 
R; 
will differ by less than A from the right-hand side of (9). Also, if 
we write L,(7') for the total length of those intervals in —T <t < T 
in which {(o+-it) belongs to R;, the numbers 
- 
Vi V{C(o+it)}{1—d(t)} dt and p3 Y'(z;) 


. 


T 


u®) (11) 


will differ by less than A for any value of 7’. Now by Theorem A 
the right-hand expression in (11) tends to (10) as 7’ > 00. If we com- 
bine this with the fact that A = A(é) > 0 as 8 > 0, we obtain (9) and 
so the proof is completed. 


THEOREM 2. Let o > 3 be fixed, and let V'(z) > 0 and (z) be con- 


tinuous functions of z such that 


im 9 = V{l(o+it)} dt< a (12) 


-T 
and O(z) = o{¥(z)} as |z| > 0. (13) 
Then the mean value (1) exists and its value is given by (5). 

Proof. We prove that (z) satisfies the conditions of Theorem 1. 
To see this we write, for an arbitrary step-function ¢4(t) which takes 
the values 0 and 1 only, 


- 
— 1 f[ 

umn oT | f(t) dt 
a 


If |O(z)| < K when |z| < M and |®(z)| < p¥(z) when |z| > M, we 


ha ve clearly T 
lim om | IO{L(o-Lit)} |4(t) dt < Ka+p0, 
a | 


* That some lines in R are left out by this division is of no importance. 
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where C denotes the left-hand side of (12). Consequently (8) will be 
satisfied if we first take M so large that the corresponding p < «/2C, 
and then, with the corresponding value of K, take » = n(e) = «/2K. 

4. Theorem 2 contains the theorem enunciated in the introduction, 
viz. that (3) implies the existence of (2) for any positive k < m, and 
it also yields for the value of (2) the expression 


G,(o) = [f \2**F(e;0) dude. 


We also obtain the existence of the mean value 
- 
L wee 7 | {(o+it)—a|** dt = C,_,(c) 
p 
and the formula 7 
kal®) = |} |z—al?*F(z; 0) dudv 
x 
for any complex a and any positive k << m. Other choices of ®(z), 
which lead to interesting mean values, are @(z) = log|z| and 
M(z) = log|z—a|, but these cases are not covered by Theorem 2, 
which deals only with continuous ®(z). It could, however, be shown 
that also in these cases the mean value (1) exists and is given by (5). 
Finally, it may be useful to observe that considerations similar to 
those of this paper can be made for other functions than the zeta- 
function. We have chosen this function because in this case the mean 
values (2) are of particular interest; the paper is rightly looked upon 
as an investigation concerning generalized almost periodic functions, 
although we have found it natural not to make any reference to the 
theory of these functions. 
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SOME SUMMATIONS OVER THE LATTICE POINTS 
OF A CIRCLE (I) 


By A. L. DIXON (Ozford) and W. L. FERRAR (Ozford) 


[Received 12 December 1933] 


1. Introduction 
1. An identity* which, with its generalizations, has attracted some 
attention in recent years is the following: 
If r(n) denotes the number of ways nf expressing n as the sum of two 
squares, X is not an integer, and v > 3, then 


MF r(n)n- IT {2m,/(nd)} = ()}"5-r(n)(A—ny-2. (1.11) 


n=0 n=0 
In the present paper we consider series that differ from (1.11) in 
having the Bessel function K, in place of J,. Starting from a known 
summation formula we are able to prove the ig: 


wii s sti K,,{2m/(an+-ab)} _ _ pha-p) ¥ of K,_,{2m/(bn+ab)} 
n=0 (n-+-b)i n=0 (n+-a)#- aad 
We then show, by the introduction, where necessary, of non-con- 
vergent series and by the use of a special process of assigning sums to 
them, that this formula contains most of the known and many new 
results concerning series whose nth term involves the number r(7). 
All previously known results about the non-convergent series met 
with in the present paper are in terms of Riesz or Cesaro sums. In 
§ 7 we prove that our process of assigning a sum to a series gives the 
Cesaro sum when such exists. 
In Part II we prove that the non-convergent series of Part I have 
Cesaro sums and, in proving this, we establish the fact that, if A is 
not an integer, ‘ 


r(n)n8J,,{27,/(na)} 


n=1 ' 
is summable (C,j) when j7 > 0 and s < 4j—}. This result differs 
from previously published results in two respects: 

(i) Previous work has been confined to series of the type (1.11), 
whereas in the statement of our theorem yp and s are independent 
numbers. 

* References to points mentioned in the introduction will be given in later 


sections. 
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(ii) Even when the work is confined to the series (1.11), the scale 
of summability has been proved only for an unduly restricted range 
of values of v; this scale we have extended to its natural limits. Our 
proof of it avoids the restriction which arises from lack of knowledge 
of the least possible value of @ in the equation 
n<N 


> r(n) = zN+O(N%). 


n=0 
2. The 7(n) summation formula 
THEeorEM A. Jf 0< a < 8, f(y) is real and of bounded variation 
in ( x, B), then ~ 
y r(n)pHf(n—0)+f(n+0)} = m 2 r(n) | feJof2nV(ny)} dy, (2.1) 


<n<p 4 
provided also that, for special values of « and B, the following modifica- 
tions of the left-hand side are made: 
(i) when n = “* the coefficient of r(B) is : be 4 f(8—9), 

(ii) when n = a + 0, the coefficient of r(x) is to be 4f( a+0), 

(iii) when n = a = 0, the coefficient of r(0) is to be f(+-0); 
and r(Q) is to be reckoned unity. 

This theorem, given in Landau’s treatise,* is the basis of our 
subsequent work. 

The formula may be written in the equivalent form 


> dr(nytf( f(n—0)+f(n+0)} 


x<nw<p 


3 
m | fly) dy +P(B)f(B)—P(a)f(a)—FB)S'(B) + Rofo) + 
B 
+ | Bw fy) dy, 
where ‘ 
Py) -"S'r(n)— Iry)—ny = yt Y r(nm-AH2ainy)}— (y > 0), 


and is O(vy) for large y;T 
msy = 
Ply) = > r(n\(y—n)—}ry* = 2 $s ny)} 


n—0 
and is O(y*) for large y. 
The modification (iii) of (2.1) when a = 0 is due to the fact that 
* Vorlesungen tiber Zahlentheorie (Leipzig, 1927), ii. 274 (Satz 559). 
+ The various refinements of the order of P(y) are not needed here. 
3695-5 E 
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P(0) = 1, while every term in the expansion of P(y) as a series of 
Bessel functions tends to zero as y tends to zero. Moreover, since 
the series for P,(y) is absolutely convergent, it is a relatively easy 
matter to obtain sufficient conditions* to justify the substitution 
B =o. We shall, in fact, apply (2.1) with’8 = oo to functions f(y), 
not necessarily real, such that 
Sy), fy), f’(y) aret bounded in (0,00) | (2.2) 
and, for large y, are Ofexp(—y*)}, where k > 0. J ie 
For such functions the proof of the ae of (2.1) with B = o 
is, either directly or by the use of the conditions referred to above, 
an easy matter. 
Assuming, then, that f(y) satisfies the conditions (2.2), we have 


bs r(n)f(n) = Sr (n) | so J(nt)} dt. (2.3) 

3. The series > r(n)n'¥K,,{27,/(nx)} 

Throughout this section R(vx) > 0. 

3.1. In our first application of the summation formula we put 
f(t) = OK, 2mi(tx)} — [R(u) > O]. 

The value of (0) is }(avx)-¥T'(). 
For large values of t, 
f(t) = Oftte-te-27vi@)} 


and the conditions (2.2) are satisfied. Hence, by (2.3), 


2 
io @] 


> r(n)n¥-K,{27)/(nx)} = 2 ¥ r( n) | HK, {2c4/(at)}Jo{2m/(nt)} dt. 
(3.11) 


n=0 n=0 0 


But 
| HK {2mrs/(at)}Jp{2rr/(nt)} dt = 2 | O+HK (27x) Jy(2arOa/n) dO, 


0 0 


* There does not appear to be any available reference to such conditions. 
On the other hand, the substitution 8 = oo in the corresponding d(n) formula | 
is readily accessible, e.g. Wilton, Quart. J. of Math. (Oxford), 3 (1932), 26-32, 
at 29. A comparison of his Theorems I and II shows that the conditions which 
govern the d(n) formula suffice for _ r(n) formula. See also Dixon and 
Ferrar, Quart. J. of Math. (Oxford), 2 (1931), 31-54, at 52-4; this proof of 
the d(n) formula is easily carried over to the r(n) formula. 

t+ At the point y = 0 the derivatives need only be right-hand derivatives. 
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which, by a known formula,* is equal to 
dar ?-Pate(n+-ax) eID (w+ 1). 
Hence, when R(y) > 0 
0 te 1 
¥ r(nynieK {2my(nze)} = Ht >} rn) __ (3.12) 


n=0 Qrht1 4 (n +a)e+ 





The proof may also be carried out by starting from the second 
member of (3.12) and using the summation formula (2.3): the neces- 
sary Bessel-function formula is given in Watson, loc. cit., 434 (2). 
3.2. We now obtain a modified form of (3.12) which will hold when 
R(u) > —1. 
When R(x) > 0, (3.12) is true and 


S r(nyn-#d = 40(ut+-1)y(u+1) = 4 SF n-# SF (—)"42n—1)21 


n=1 n=1 n=1 
Hence, we have 


D(u+1) Tu) 


2 ¥ r(n)ni#¥K,{2m,/(nx)} = anal ae 


n=1 


= 


To = 1 
47 > OY ae Seno —sa}+ £e+1)nu-+1)]. 


\(c-+-n)jet 
(3.21) 
This result is, by analytic continuation, true when R(u) > —1. 
When p = 0 we obtaint 


‘ : l 
2S r(n)Kg{2eq(ne)} = — += a rin) 7} + 


+im| — Tu) delet 10-1) meth. (3.22) 


p—>0 art bot qre tt 


After using the functional equation 
mT(1—s)C(1—s)n(1—s) = m'*1(s)e(8)n(s) (3.23) 


as a first step in evaluating the limit required in (3.22), we obtain 


r(n) Ko{27,/(nx)}—log(*x) — ash 2¢’(0) + 2n'(0) 


ats Dt a a aia 


G. N. Watson, Theory of Bessel Functions (Cambridge, 1922), 410 (2). 
+ Watson, Quart. J. of Math. (Oxford), 2 (1931), 298-309, at 300, sets out 
a similar problem in full detail. 





52 A. L. DIXON AND W. L. FERRAR 

The above process of extending the region of validity may be con- 
tinued ad libitum. The second extension gives, as a special case, a 
formula of some interest. We write the last line of (3.21) as 


vel (n+l io, { 1 1 x(u-+1)) 
gt | >) at nh+2 jv 
+40(u+1)9(u+1)— 40 (H+ 1)C(u+2)n(u+2)], 
and the formula so obtained is valid when R(x) > —2. When we 
put » = —1, we get* 


i) 


2at ¥ r(n)n* Ky {27,)(na)} 
[z+ log 7 | —logx-+2¢'(0)—29/(0)+ 
+-ae{—log(n®x)-+1—2y—2£'(0)-+2n/(0)}. (8.25) 
As one might expect, (3.24) is merely the result of differentiating 
(3.25) with respect to 2. 
4. The series § r(n)(n-++-b)-*“K,,{27,/(an+ab)} 
Throughout this section R(va) > 0,b > 0. 
4.1. When b > 0, the function 
S(t) = (+6) K, (27, (at-+-ab)} (4.11) 
satisfies, irrespective of the value of yu, the conditions (2.2). Accord- 


ingly, by (2.3), 


; 
2 rn )(n-+-b) 4#K, {224 (an+-ab)} 
n 


2) 
oo . 
= 7 > r(n) | (¢-+6)-#K,{27y(at+ab)}Jo{27r/(nt)} dt. (4.12) 
n=0 0 
The integrals in (4.12) are of Sonine’s type.t We write a = a2, 
b = B*, t = & and so obtain as the coefficient of r(n) in the second 
member of (4.12) 
00 
‘ 21 Q2)-twK fs (92. R2 ‘ 
2a | (0°?) *#K, {2ara/(6°+- B?)}Jg(27rOvn) 0 dd 
0 
— Rl-Hy-H(y2-1)-14+eEK. f5 M2 
= B)-Ha-#(a?-+-n)-? 4K, ui2mBy (x?-+-)}. 
* The odd terms at the end require, of course, considerable arithmetical 
detail for their evaluation. 
ft The formula we use to evaluate the integrals is given by Watson, loc. cit., 
416 (2). 
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Hence we may write (4.12) in the form 


aie ¥ r( n)(n-+-b)-#K, {27,/(an+-ab)} 


n=0 
= bi0-») > r(n)(n-+a)-10-K,_,{2/(bn-+-ab)}. (4.13) 
n 
4.2. The formula (3.12) is, essentially, a particular case of the 
formula (4.13), though we find it convenient to give independent 
proofs. If, in (4.13), we suppose R(x) < 0 and take the limit, term by 
term, as b + 0, we obtain (3.12) save for a trivial change from yp to —p. 
4.3. A special case of the formula (4.13) was discovered by Rama- 
nujan. It is given, with an indication of how it was obtained, at the 
end of a paper by Hardy.* 
In (4.13) put » = }. It becomes, since K,(z) = e~*,/(7/2z), 


e-27v(bn+ab) — a a aaa (4.31) 
n 


which is precisely Ramanujan’s formula. 


5. The formula p > r(nje—"*?* — p- ¥ r( r(n)e—"7!P* 


n=0 n=0 
This formula, which may be obtained, when R(p?) > 0, by puttingT 
S(t) = exp(—p*at) 
in (2.3), is also obviously the result of squaring each side of Poisson’s 


formula 


r. > en ‘7?’ — 9- $ >) en *a/p* 


n2— 0 

The formula is ied aaa’ with the bulk of our present investi- 
gations: we note it here (a) for the sake of completeness, (b) because 
the d(n) identity which corresponds to it has certain points of interest 
that we hope to consider later. 


6. Deductions from the K series when 6 > 0 
Throughout this section we assume b > 0 save in the definition of 
§ 6.2, where b > 0. 
1. The formula (4.13) has been proved subject to the hypotheses 
- 0, R(va) > 0. We now examine what happens when R(va) = 0. 
[The details are somewhat lengthy, but the reader may see for 
himself how the actual formulae run if, at this point, he puts va = ivA 


* Quart. J. of Math. 46 (1915), 263-83. 
+ Watson, loc. cit., 393 (i), gives the values of the resulting integrals. 
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in (4.13), notes that K,(ix) is a multiple of J,(x)—iY,(x), and 
assumes that the series which arise have, for the most part, Cesaro 


sums. | 
n (4.13) put va = a+ivA. Then,* for « > 0, 


(a-+-avaA)e ye ¥ r(n)\ n)(n-+-b)-4#K , {207(a-++-iva),/(n+6)} 
n= 
= 510-») > (n){n+(a-Fivd)?}U YK, _s{2arvb,/(n+ 02 —A-+ 2aiva)}. 
nt 
(6.11) 
From the asymptotic formula 
yr fo) — —Zy—4 
K,,(2) = O(e*z>), 
we see that both series in (6.11) remain absolutely convergent when 
a = 0, provided that » > 3. For such values of » and for any fixed 
positive value of A the series converge uniformly in a > 0, and we 
have, by making « — 0 in (6.11), 
elpmidin 5 r(n)(m+-b) 4K, {2ari,/(An-+-Ab)} 


n=0 


n<r 
= b10-weime—-Di } r(n)(A—n)*e- g i,/(Ab—nb)}+ 
n=0 
+bi-») ¥ r(n)(n—A)Me- YK, -s{2m,/(nb—Ab)}. (6.12) 
n>A 


On using the formula 
et K (ai) = —hriH? (2), (6.13) 


we see that (6.12) may be written as 


Ate >: r(n)(n-+-b)-# HP ){2ar,/(An-+-Ab)} 


n 


a = pia- ) ae A—n) jie a. i{227,/ (Ab—nb)}+ 
n=0 
+ (2i/7)bi0-») > )(n—A)#e- Om... a4/(nb—Xb)}. (6.14) 


n>A 


The real parts of (6.14) give 
die s r(n)(m+-b) J, {2274/(An-++-Ab)} 


n=0 


— pia-p) s r(n)(A—n)KH-DJ,_s{27,/(Ab—nb)}. (6.15) 
n=0 


* We change from K,_, to K,_, and so avoid considerable arithmetical 
labour when, as later, we put a = 0. It is the functions J, Y,—,; and not 


w— +u— 
J,_,» Y,_, which appear in the sequel. 

















SUMMATIONS OVER LATTICE POINTS OF A CIRCLE 55 
Before proceeding further it is of some interest to examine the rela- 
tions between the last formula, (6.15), and the well-known identity 
(1.11). The formula (1.11) is merely the result of putting 6 = 0 in 
(6.15). In fact, when we expand both sides of (6.15) as power series 
in 6 {the first side by Lommel’s expansion* and the second by the 
expression for J,(x) as a power series in 2], the fact that the two 
coefficients of 6” are equal is expressed by (1.11) not only when 
m = 0 but also when m = 1, 2, 3..... 
The imaginary parts of (6.14) give 


du S r(n)(n-+-b)-MY, {2n/(An-+2b)} 


n=0 


— pia-») yer r(n)(A—n)#-DY_ _,{27r,/(Ab—nb)}— 


n=0 
— (2/ar)bt@ -») ¥ r( )(n—A)H-D Kf 22r4/(nb—Ab)}. (6.16) 
n>dA 
6.2. So far we have considered (6.11) with « = 0 only when p > 3. 
But, irrespective of the value of yu, the second series in (6.11) is, in 
general,} convergent when « = 0 and its sum is continuous in « > 0. 
When «a > 0, the first series in (6.11) is convergent; when a = 0, it 
is not necessarily convergent. We may, however, consider that this 
first series has, for « = 0, a conventional sum equal to the limit of its 
sum as «+ 0, i.e. save for certain exceptions, the convergent sum 
of the second series in (6.11) when a = 0. 
We now give a formal definition of this conventional sum, but we 
defer the proof of its consistency with ordinary convergent sums and 
with Cesaro sums until the next section (§ 7). 


DEFINITION. Given a series 
> a, K,,{(a+%B)y(n+b)} (6 > 0) (6.21) 


which converges for « > 0 and defines a function F(x); given also that 
F(a) >l as «—> 0; then I is said to be the Abel-Bessel, or AB, sum of 


the series ya, K, {iB (n-+b)}. (6.22) 
A series $a,J,{By(n+b)} is said to be summable AB when the 

series > a, H\?4B,/(n+-6)}, that is 
(2i/ar)etu7’ ¥ a, K,{t8,/(n+5)}, (6.23) 


* Watson, loc. cit., 140. 
+ The exceptions are given at the end of the present sub-section. 
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is summable AB: the sum of the J series is the real part of the sum 
of (6.23). 

There is a similar definition for Y series. 

Consider now the equation (6.14), a form taken by (6.11) if « = 0. 
When the left-hand side is interpreted as an AB sum, the equation 
(6.14) is true* 

(i) for all values of p, if r(A) = 0; 

(ii) for values of » such that R(u) > 1, if r(A) + 0. 

In such cases the right-hand side of (6.11) is continuous in a > 0 and 
the limit of its sum as a 0 is obtained by putting « = 0 in the 
separate terms of the series. 

On the other hand, if r(A) 4 0 and R(u) < 1, the right-hand side 
of (6.11) is not amenable to the limiting process required in the 
definition of an AB sum. Actually, this right-hand side contains 
a term (n = A) which 

(i) when » = 1 or R(u) < 1, tends to infinity as « tends to zero; 

(ii) when » = 1++7y and y + 0, has no limit as « tends to zero; 
for if we put a?+2aivA = ¢ and consider the expansion near t = 0 
of H4-DK,,_,{27,/(bt)}, the dominant term is }(¢7vb)+—"T(1—y) for 
p ~ 1 and —log{z,/(bt)} for p = 1. 

In either case the left-hand side of (6.14) cannot have a finite AB 
sum, though it is interesting to notice that when » = | the real part 
remains finite while the imaginary part does not. 


7. The consistency of 4’ sums with Cesaro sums 
Throughout this section b > 0. 


7.1. We now prove the following theorem. 


THEOREM B. Jf b > 0 and the series 
> 4, K,,{i8,(n+b)} (7.11) 


is summable (C,k) where k > 0, then its sum is given by 


lim > a, K,,{(~+%B)/(n+6)}; (7.12) 
a0 


that is, if (7.11) is summable (C,k) where k > 0, then it is also sum- 
mable AB and the two sums are identical. 


* (A) 0 when (i) A is not an integer (or zero), (ii) when A is an integer 
that cannot be expressed as the sum of two integer squares: 7(A) 4 0 when 
A is zero or is an integer which can be expressed as the sum of two integer 
squares. 
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The proof is simple but somewhat long. There are two essential 
points. The first is the asymptotic expansion 
K,(2) ~ (F)Pe b ome (7.13) 
where (u,m) = D(u+m+$4)/P(m+ 10 (u—m-+ 4). 
The expansion is valid when |argz| < 37. The second point, which 
selects the leading term of the asymptotic expansion and establishes 
an Abel theorem in connexion with it, is contained in the following 
lemma. 
LemMaA 1. If > a, ts summable (C,k) where k > 0, then its sum is 


given by lim y a, e-¥vinsd) 
. } 
y>0 


m=0 


Since the (C,k) sum of a series is not altered when k is increased, 
it is sufficient to prove the theorem when k& is an integer. Actually, 
we prove it here for k = 1; the details when & is an integer greater 
than unity are a little more elaborate while the process is much the 
same. 

Let 8, = dg+a,+...+4,, o, = 8+8,+...+8,, and suppose that 
> a, has a (C,1) sum s. Then 

o, = ns+o(n), 8, =o(n), a, = 0(n). 


Now, by a well-known formula in integral calculus, 


i 2) 
e-uvinst) — 2 [ ¢-a-arinsyte* de, (7.14) 
Vir | 
0 
and, by partial summation, 
N 9 N m - 
> a, yv(n+b) — —- > 8, [ e-v*—y"(n +b)/dar*(] — e-v'lA4zx*) ‘ie ra 


[ O(N eV +dx* dae, 
the last term being O{Ne-vv\V+)}, . 


Thus, 


A 
VS a. e~yvay) — 
la TR 


2 © 
$| v9 


n=0 


_ oO 
—7—4% b)/4a2 _ytlay? ss 
> fe x y(n-+b)/4x*(] —¢ y*/4a ) dx 
n=0 0 


i | e-2*—vinsyac"(] —e-v'l8e")2 dr, (7,15) 
) 


0 


4) 
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For convenience, we write 


dy, = e-o*—'(n-+b)/4a*( J —e-v' lay, 


and note, as a preliminary, that 
a ~ 
= > 0 | by da = er’, 7.16 
Nar 2 $n ( 
n=0 0 


Given any positive «, we can find an m, such that, for n > mp, 


(o,,/n) = s+te,, where |e,| <e, 
and so we can find a constant K such that, for all n, 
(c,/n) = 8+K,, where |K,| < K. 
We have, therefore, from (7.15), 


io 2) @ 


© ie) 
oo «© ° No - © ° 
3Vrr > a, e-Uvint) = 8 > n | d, dx + > nK,, | d, dx + ¥ ne, | py, dx. 
n=0 n=0 n=0 n>Ne : 
0 0 0 
Hence, on using (7.16), 


i 6) 
i. 


n=0 n=0 


(oe) 
K No . wo 
dvr} > a, e-UVrt+)—ge-vv) < Kn | ¢, dx +e> n | d,, dx. 
: n=0 
0 0 


That is, 


No 


2) 
eo 
Zz ay e-vvins)_g| < |s(1—e-¥v) |+-ee-¥v0+ (2K, vi) > n | Pp dx ? 


n=0 n=0 
0 
and, , being fixed, the last term tends to zero with y, as may be 
seen at once on using (7.14). The truth of the lemma is now evident. 
7.2. We next establish a lemma which, in due course, enables us 
to use Lemma | in the proof of Theorem B. 
Lemma 2. If k is a positive integer and 
D> 4, K, {iB (n+6)} (7.21) 
is summable (C,k), then so is 
> a, (n+b)-* exp{—iB,/(n+b)}. (7.22) 
As before, we prove the lemma for k = 1 as a sufficient indication 
of the proof when k is any positive integer. Write 
C, = a, K, {iB (n+b)}, 
In = (n+b)* exp{—iB,/(n+5)}/K,{iB,(n+b)}, 
and suppose > c,, has a (C,1) sum. 
From the asymptotic expansion (7.13), we have 


fn = Ag +Ax(n+b)+*+...+Ay(n+b)*+ O(n-+b)-2, 
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where the A are independent of n. Moreover, since > c,, has a (C, 1) 
sum, ¢,, is at most O(n). 

Hence* > c,f, may be considered as the sum of five series, 
> A,,¢,(n+b)-™ (m = 0, 1,..., 4) each with a (C,1) sum, together 
with a sixth series whose terms are O(n-). It follows that > c, f,, is 
summable (C, 1), i.e. (7.22) is summable (C, 1). 

7.3. Proof of Theorem B when k = 1. 


We first define an auxiliary function ¢,(z) by means of the 
equation } $.( ) 
_ wa Se p,m 
$,(z) = K,,(z) (=) e* >, (22) (7.31) 
If R(z) > 0, $,(z) = O(|z|-*), as may be seen from (7.13). 
In or a to Theorem B when k = 1 it is sufficient to prove 
() that Lang fiByin+0}— (C,1) 
= lim > a, ¢,{(«+%B),/(n+)}, 
a->0 
and (ii) that 


exp{—iB,/(n+5)} 
Ss ay (2Bi)*+™(n-+-b)t+im (C, 1) 


: b)} 
=tim > a, exp{—(a+7B),/(m+ 


*(2(a-+ spy} (m+ b) Ham 
when m = 0, I 
In the first place, (i) follows from elementary considerations. By the 
hypothesis that (7.11) is summable (C, 1), a,, is at most O(n*/*), and we 
see from (7.13) that the nth term of each series involved is at most 
O(n-*), and so the second series is uniformly convergent in « > 0. 
In the second place, Lemma 2 shows that 


< exp{—iB,/(n+b)} h 
ia m+ 2, (m = 0, 1,..., 4) 


(n+b)im+t 

is summable (C,1). By Lemma 1, its (C, 1) sum is given by 

um 50, sri lo tele 

a m+4 aa ” (n+b)im+t 
—(a-+ iB). 
= lim ¥ a, SPO +PWln tO) 

a0 & {2(x-+-aB)}"+#(n-+-b)im+ 

and so (ii) is proved. 
This proves Theorem B when k = 1. The proof when k is any 

positive integer may be written out in much the same way. 


* If > c, is summable (C, 1), then so is > c,(n+6)~* for any positive s. 
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7.4. We shall return later (Part I1) to the problem ‘When are the 
non-convergent series in (6.12) [the equation which led to the intro- 
duction of summability AB] summable by Cesaro means?’ For the 
moment we note that the exceptional cases given at the end of § 6.2 
show that, when r(A) ~ 0 and R(u) < 1, the series 

> 1(n)(m+b) + K, {277i (An +-Ab)} 
is not summable by Cesaro means. If it were so summable, its sum, 
in virtue of Theorem B, would be evaluated from the other member 
of (6.12) and that contains an infinite or an undefined term. 

Oppenheim* and Wilton} have each noted cases where similar 
series were not summable by Cesaro means. 

7.5. We conclude this section by stating a theorem which sum- 
marizes most of the positive results (as distinct from the negative 
results noted in 7.4) of §§ 6 and 7 in so far as they concern the r(n) 
type of series we are considering. 

THEOREM C. Jf b > 0 and the series 

AMF r(n)(n-+-b) HH 242mr/(An-+Ab)} 
n=0 


is summable by Cesdro means, then its sum is given by 


ns A 
biA-/#) ¥ r(n)(A—n)U DA {221/(Ab—nb)}+- 


n=0 


+ (2%/2)btG-) s r(n)(n—A)MH-VK, _{ 2ar,/(nb—Db)}. 
n>dX 
The real and imaginary parts of these expressions are noted in the 
formulae (6.15) and (6.16) respectively. 
The corresponding results when b = 0 are considered in § 8. 


8. Deductions from the K series when b = 0 
. There are a number of formulae, mostly new, which follow 
of om < results of § 3. The most interesting aret 


. - / 
(i) > r(n)nied_,{2my(nd)} (AB) 
n=0 
’ <A 
MYT (+1). " 
— _ sin pur 
ate (A—n)e+1 
n=0 
* Proc. London Math. Soc. (2) 26 (1927), 298; Theorem 4. 
+ Mess. of Math. 58 (1928), 76. 
t The notation (AB) implies that the sum of the series is an (AB) sum: 
if r(A) ~ 0 (i) and (ii) may not be true, but they then indicate the nature of - 
the divergence of the first series in each of these equations. 


a (8.11) 


> 
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(ii) Sr(n)n“J{2ny(nd)} (AB) 


n 


_ AD (w+) r(n) 
sin pa > 


mete \u+ 


(n—A)jeH 
when R(u) > 9; 


(iii) ¥ r(n)[¥of2m(nd)}+(2/m)Ko(2m(nd)}] (AB) 


= 7 ‘dr Ncats = (8.13) 


n 
8.2. We have proved the Simin (3.12) on the hypotheses 
R(u) > 0 and R(vxz) >0. Suppose now that R(z) > 0, that 
Vx = a+ivd, and that a tends to zero. We get 


io 2 


 Sr(nyninK (2mrif(nd)} (AB) 
n<A —_lumt = mi 
e et [= +S ile | (8.21) 


2a +1 A—n)e+1 n—A)e+t 
0 A 


the second member being convergent save when r(A) + 0. 
In this we put 
K,,{2mi,|(nd)} = — Jarre tot {2rr,/(nd)} 


and so, when R(u) > 0 and r(A) = 0, we obtain (8.12) together with 


5 r(n)ni¥Y {2zr,/(nA)} (AB) 


n=0 


— MED (w+) RA v(m) =< r(n) » 
= = > Gaara omer > Gal (8.22) 


n=0 
When we combine these results by means of the formula 
J_,(z) = J,(z)cos yr —Y,,(z)sin pr 
we obtain (8.11). 
We have thus proved, as may be seen by putting » = —v in (8.11), 


NYS (nn, {2m/(nd)} (AB) 
n=0 
= {9’/P(v)} WS r A—n) (8.23) 


provided that R(v) < 0 and r(A) = 0. 
This identity, for values of v other than those we have yet con- 
sidered, is well known and has been frequently discussed.* When 
* Wilton, ‘The lattice points of a circle; an historical account of the 


problem’: Messenger of Math. 58 (1928), 67—80, gives an excellent introduction 
to the topics connected with this identity. 
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v > $ the first series in (8.23) is convergent and its sum is given by 
the second series.* When 4 < _v < } Wiltonf has proved the formula 
(8.23) save that he uses summability (C, -+-v-+8) instead of the (A B) 
method of summation. 

We conclude this sub-section by indicating how (8.23) is proved 
when 0 < v < 3. Instead of using (3.12), which requires R(u) > 0 
we use its extension, namely (3.21), which is valid when R(x) > — 
In (3.21) put » = —v and we get, on noting that K_,(x) = K,(z2), 
25 r(n)n— K,{2n)(na)} = —— ch tl wl 


n=1 


_ ee 
vl hv oi 


+D(1—v)2ta-»| y r(n){(w-+-n)*—n?-3} + 4¢(1—v)(1—v)]; (8.24) 
n=1 

this is valid when R(v) < 1 and R(vx) > 0. On writing, as before, 
Vx = a+ivA and then making « tend to zero, we have, as one of the 
two results obtained by taking real and imaginary parts, the equa- 
tion (8.23). 

We have thus proved (8.23) when r(A) = 0 and R(v) < 4. If we 
anticipate a result of Part II, namely, that the first series in (8.23) is 
summable (C, k) for k > 4—v, we may, in virtue of Theorem B, write 
our result in the form 


@ 


YY r(n)n- ST {27,)(na)} (C,4—v+e) 


n=0 
n<A 
= {r’/T(v)} ¥ r(n)(A—nyr 
n=0 
whenever 7(A) = 0 and v < }. 
9. Combination of series which have AB sums so as to 
give convergent series 
In many cases series with AB sums can be so combined that the 
resulting series is convergent when the terms are properly grouped. 
For example, if we take (6.15) for b > 0, and also for b = 0, and 
subtract one result from the other, we get 
AHS r(n)[(m+-b)- ed, {2ar/(An +-Ab)}—n- 1S {2ar,/(An}] 


n=0 


n<r 
= bi0-») Y r(n)(A—n)he-YT, > }— 


n=0 


—for-4/ uy" )(A—n)e- 


* There are several known ways of proving this: our own method is indicated 
in a former paper, Quart. J. of Math. (Oxford), 2 (1931), 31-54; at 33. 
+ Proc. London Math. Soc. (2) 29 (1929), 168-88, Theorem D. 
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and now the series on the left is convergent when p > —}, as is 
easily seen by carrying on the expansion of (n-+-b)~#J,{27,/(An-+-Ab} 
to two more terms in ascending powers of b. The AB sums of § 6 may 
be treated in the same way. Further, if we notice that in the defini- 
tion (6.2) of an AB sum, B may be zero or a complex number, we see 
that we can derive (3.12) and (3.21) from (4.13) by putting a = 0 
and then applying our present process, 











“ A NEWTONIAN EXPANDING UNIVERSE 
By E. A. MILNE (Ozford) 
[Received 7 March 1934] 


1. The phenomenon of the expansion of the universe has usually 
been discussed by students of relativity by means of the concept of 
‘expanding space’. This concept, though mathematically significant, 
has by itself no physical content; it is merely the choice of a particular 
mathematical apparatus for describing and analysing phenomena. 
An alternative procedure is to choose a static space, as in ordinary 
physics, and analyse the expansion-phenomenon as actual motions 
in this space. Moving particles in a static space will give the same 
observable phenomena as stationary particles in ‘expanding’ space. 
In each case the space is a construct built up by the mathematician 
out of observations that could in principle be made; it is built up 
around the matter in motion according to certain rules. The formu- 
lation of the relevant laws of nature depends on the rules adopted, 
and the laws will be quite different if different rules are adopted, as 
I have elsewhere* explained. The alternative procedures have been 
tersely described in a recent paper by 8. R. Milner.t He explained 
that we can either modify our geometry in order to retain 6 { ds = 0 
as the paths of free particles, or retain Euclidean geometry and Min- 
kowski space-time and modify the variational principle by weighting 
the elements of path ds with appropriate invariant weighting factors. 
Einstein’s general relativity adopts the first procedure; in my recent 
treatment of the cosmological problem I adopted the second pro- 
cedure. In neither case has the space constructed any physical signifi- 
cance—no content attaches to the phrase ‘the space of nature’—but 
the second procedure has the advantage that it employs the space 
commonly used in physics. The equations of motion I obtained for free 
particles moving in the presence of a certain distribution of matter in 
motion were derived by the straightforward method of making them 
invariant under the relevant transformations, but, in confirmation of 
Milner’s exposition, Dr. A. G. Walker has recently expressed themt 

* Zeits. fiir Astrophys. 6 (1933), 29. 

+ Proc. Roy. Soc. A, 139 (1933), 349 (§ 2). 

t In course of preparation for publication. The acceleration formulae in 
question were given in Zeits. fiir Astrophys. 6 (1933), 56, equations (70). 
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in the form of a principle of least action, and obtained the corre- 
sponding weighting factors. 

Milner has remarked that formulae obtained by the one method 
should be convertible into formulae obtained by the other. For 
example, the formulae describing the velocity- and density-laws for 
an expanding universe obtained by the methods of the general theory 
of relativity should be equivalent to similar formulae obtained when 
the matter-in-motion is described in flat, static (i.e. Euclidean) space, 
on a suitable re-formulation of the relevant laws of nature. 

In this paper I show how the same locally observable results can 
be obtained from elementary Newtonian theory (using flat, static 
space, Newtonian time, and the Newtonian dynamics and law of 
gravitation) as are given by Einstein and de Sitter’s well-known rela- 
tivistic model* of a universe in flat, expanding space and the relativistic 
theory of gravitation. It will be shown that the latter model corre- 
sponds to a Newtonian universe in which every particle has the 
parabolic velocity of escape from the matter ‘inside it’ as judged by 
any arbitrary observer situated on any particle of the system. It will 
further be shown that the governing equations in the Einstein-de 
Sitter case, namely those describing the behaviour of a particle with 
fixed coordinates, are identical in form with the Newtonian equations 
governing the distance of a particle from the observer as functions of 
the time. The identity is exact—no approximation is involved, nor is 
any neglect made of inverse powers of the velocity of light. These 
results will be extended to a general class of relativistic universes 
and the corresponding Newtonian universes in a joint paper by the 
author and Dr. W. H. McCrea. 

In the Newtonian cases the symbol ¢ occuring in the differential 
equations and their integrals denotes Newtonian time. In the rela- 
tivistic cases it denotes ‘cosmic time’, i.e. the time kept by a clock 
moving with the particle concerned. In the Newtonian case such a 
clock keeps the same time as the observer’s clock, assuming the usual 
definition of simultaneity by means of light-signals; in the relativistic 
cases, it can be shown that the ‘cosmic time’ of an event does not 
coincide with the epoch assigned to it by a distant observer, using the 
same definition of simultaneity. Thus the two identical sets of 
differential equations have different interpretations in the two cases. 
This shows that the contents of the two theories are not given simply 

* Proc. Nat. Acad. Sci. 18 (1932), 213. 
F 
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by the resulting differential equations, a fact which seems of some 
significance in view of the opinions of some thinkers that it is the 
differential equations alone which convey what we are entitled to 
believe about nature. Apart from this question of interpretation, the 
relativistic and Newtonian theories as regards models of the universe 
are indistinguishable in their predictions of local phenomena—and it 
is to be remembered that hitherto the only phenomena astronomically 
observed must be counted as strictly local phenomena, for no recession 
velocities exceeding say one-tenth of the velocity of light have yet 
been measured, and up to this limit there is little difference between 


the predicted phenomena. 


2. It seems to have escaped previous notice that whereas the theory 
of the expanding universe is generally held to be one of the fruits of 
the theory of relativity, actually all the at-present-observable 
phenomena could have been predicted by the founders of mathe- 
matical hydrodynamics in the eighteenth century, or even by Newton 
himself. The velocity of light, c, does not enter into the formula 
determining the law of expansion or the relation between the rate of 
expansion and the local mean-density. This point is obscured in 
treatments which take the velocity of light as unity; actually c cancels 
out, and a knowledge of the numerical value of c is not required. All 
that is necessary is the Newtonian theory of dynamics and gravita- 
tion, combined with the hydrodynamical equation of continuity. 
Both ‘special’ and ‘general’ relativity are unnecessary; Newtonian 
relativity suffices. The resulting universe satisfies Einstein’s cosmo- 
logical principle* that all particles in it are equivalent, both in their 
local relationships and in their distant relationships, so far as these 
are unambiguously fixed by Newtonian theory. Imperfect as the New- 
tonian universe is as not being in accordance with Einstein’s principle 
of relativity, its study is valuable as affording insight into the whole 
subject. As faras present astronomical observations go it is completely 
adequate, and the actual phenomena of the expanding universe, as 
locally observed, provide no evidence for or against the principle of 
relativity in Einstein’s form, though they do provide evidence in favour 
of the Newtonian principle of relativity in its own (limited) domain. 
Lastly, study of the Newtonian universe illustrates the point that to 
use ‘expanding space’ is not the only way of dealing with the subject. 


* To be distinguished from Einstein’s principle of relativity. 











A NEWTONIAN EXPANDING UNIVERSE 67 

A Newtonian universe has been outlined by W. R. Mason,* but he 

introduced an isothermal distribution of matter and an equation of 
state. The present treatment is entirely different. 


e 

3. Velocity- and density-laws. Let us adopt Euclidian space 
and Newtonian time for all observers, and the Newtonian formula- 
tion of dynamics and gravitation. Consider a swarm or cloud of 
freely moving particles in this space; in the system we shall construct 
collisions do not occur. The problem is to find a cloud of particles, 
possibly in motion (i.e. to determine its motion and density-be- 
haviour) such that it is described in the same way as viewed from any 
particle of the system as place of observation. Einstein showed long 
agot that a static universe of this kind led to contradictions within 
the walls of Newtonian theory. We therefore investigate the possi- 
bility of constructing a homogeneous universe in which the density 
p at any point changes with the time. Since here there is an absolute 
simultaneity, there is no difficulty in defining homogeneity, and we 
therefore have p = p(t), a function of time ¢ only. Of the possible 
motions, a particular case will be that in which the direction of 
motion is strictly radial as seen by an assigned observer. (We exclude 
the possibility of rotation.) We now investigate the form of the 
function p(t) and the law of dependence of velocity on position and 
epoch. 

Let v be the outward velocity of a particle at time ¢, at distance 
r from the observer, relative to the particle on which the observers 
situated. Let M(r) be the mass contained in the sphere of radius r. 
Consider the particular case{ in which the distant particle has the 
parabolic velocity of escape from the mass contained in the sphere of 
radius r. The observer considers the material outside this sphere as 
having no influence on the motions inside it, in accordance with 
Newtonian gravitational theory; the observer, in fact, supposes that 
conditions ‘at infinity’ are compatible with this assumption. Then 
he writes down 
GM (r) 


— ae (1) 
: 


to 


vy 


* Phil. Mag. (7) 14 (1932), 386. 

+ Sitz. d. Preussischen Akad. (1917), 150. 

t This will be generalized later, in the joint paper which follows, but the 
particular case offers a better introduction and in any case requires separate 
treatment in the resulting integrations. 
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Since the mass M(r) remains constant ‘following the motion’, the 
particle will always possess the velocity given by (1) if it once 
possesses it. In writing down (1) we are not using the notion of a 
gravitational potential, here inapplicable, but are employing (1) 
simply as an integral of the equation of motion with a particular 
value of the constant of integration. Equation (1) gives 


87rG ™ 


“2. 
; 


p2 


The motion must be such that the hydrodynamical equation of 
continuity* is satisfied. This, in polar coordinates, runs in Eulerian 


notation ‘ 1 ¢ 
+5 = (r°pv) = 0, (3) 
ar’ or 
where v is a function of r and ¢ given by (2). Inserting this, we have 
a, 87G\! 
p-iFf43(-2") —0. 
dt 3 
Integrating, we have 
— 2p-}+ (247G)!t = const., 
or, choosing a suitable origin of ¢, 
—2p-!+ (247G)'4t = 0. 
This gives p = 1/6rG?. 
- . ; 2r 
Equation (2) then givest 37” 
We now verify that this is a solution of the problem. The accelera- 
tion of the particle is 


Dv _D 2r _2fv ¢ a 2r 
Dt ~=—~Dt\3t 3\t #@) ~~ oe’ 


and this is precisely the Newtonian acceleration, since 
Y » 

_ GM = _ SerGiy? _ 1 "= _ 

' 6rGt? r? 9t? 

Lastly (5) satisfies the Newtonian principle of relativity and 
Einstein’s cosmological principle. If we transform our origin to 
another of the moving particles, at distance R, where the velocity 
V is §R/t, the Newtonian formulae of transformation are 


r’=r-R, v’ = v—V, 


* Circa 1750. {+ The minus sign is also permissible, giving v = —3r 
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_ ar. ee’ 
Ve eS —_ = ——,”" 7 
3t st 3t (7) 
The acceleration also obeys the Newtonian transformation law. 


whence 


4. Discussion. Equations (4) and (5) provide a solution to our 
problem. By (5), v obeys a velocity-distance proportionality at any 
one epoch, and so this Newtonian universe obeys Hubble’s law of 
nebular velocities. If we put this law in the form v = ar, where « is 
observed, then « = 2/3¢, and by (4) 

2 


a 
p= oy (8) 
The data from the nebular velocities and distances then give for ¢ a 
value of about 1-3 10° years, and a density p of about 5x 10-%8 
gram cm.~* These are of the usual orders of magnitude given by the 
‘expanding space’ theories, as well as by the kinematic theory. In 
the kinematic theory, as I have shown,* p ~ 1/{7G#? and v = r/t, so 
that here p ~ a*/{7G. Thus the local value of the ‘age’ of the 
universe is on the Newtonian theory two-thirds that on the kinematic 
theory, and the density about one-half that on the kinematic theory. 
Present estimates of the actual mean local density of the universe 
cannot discriminate between the two. 


5. Comparison with the Einstein-de Sitter universe. We 
have seen that the Newtonian universe constructed above is defined 
by the equations 


dv? — 


GM(r) Dv _  GM(r) 
r Dt 38 PF 


; (9) 


where M(r) = $xpr’. (10) 
Equations (9) contain within themselves the equation of continuity, 
for on differentiating the first of (9) and using the second we have 
at once 
D 
— M(r) = 0. 
pi”? 
Put r = fR(), (11) 
where f is a constant particularizing the particle considered and R is 
a universal function of ¢ only. Actually, integration of 
_dr_ 2 
dt 3t 


* M.N.R.A.S..94 (1933), 7. 
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gives at once r = fti, so that R(t) = t’, but we do not need this. 
Then v = fdR/dt, Dv/Dt = f d*R/di?. Introducing these in (9) and 

using (10) we see that f divides out and we get 

2 Y 
(ji z) = —? (12) 

2d*R 82G 

Rd®” 3 


__ (Lary 
— aa 


Introduce Einstein’s constant « defined by 
Kk = 82G/c?, 
and write c dt = dr. 
Then (12) and (13) become 
| ae 
( 7) we x, 
2 Tat; 7) aii 


R dz ¢ 


Rdr 


But these are identical with the relativistic equations for an expand- 
ing universe of zero curvature with pressure p = 0 and cosmical 


constant A = 0, as given* by Einstein and de Sitter. 

Conversely, from the relativistic equations (14) and (15) we can 
infer equations identical in form with the Newtonian equations. The 
equations (14) and (15) are derived from a metric 

ds? = dr*— R?(dxz?+-dy?+dz?). 
In this space a particle is assigned fixed ‘coordinates’ x, y, z, and the 
‘distance’ r of such a particle is given by 
r= fR, 
where f is constant for the particle, depending on the particle chosen. 


Then ldr 1dR 
rdr Rdr 

Introducing these in (14) and (15) and returning to ¢, we see that c 

cancels out, and we are left with 

1 dr\? 87rG 

on 16 

( i 3° (16) 


* Einstein and de Sitter, loc. cit., give (14); (15) is given by de Sitter, 
Univ. of California Pub. Math. 2 (1933), 161. 
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2 d*r 1 dr\? 
rd@ ( 3) e 
Define m(r) by m(r) = {nrp. 


3 


D D D{ 3 _(dr\2 
Th = a= fo de on fe re 
si pi) = sr Hylrr) lear) r 


which by (17) is zero. Hence m(r) is constant following the motion, 
and (16) and (17) may then be written 

se Gm/(r) Dv ca Gm(r) 

Sie Dt re” 

which are the Newtonian equations. It follows that the two equations 
defining the behaviour of R and p in the Einstein-de Sitter universe 
are equivalent to an equation of motion and an equation of con- 
tinuity. Since they were originally obtained from Einstein’s field 
equations via the Riemann-Christoffel tensor, we have an interesting 
example of the correspondence of Einstein’s field equations with 
Newtonian dynamics and gravitation. The density p in the Einstein- 
de Sitter universe now comes out,* as in the Newtonian case, as 
p = 1/67G#?, and the ‘velocity-law’ as v = 3r/t. 

Since the time ¢ or 7 in the relativistic case coincides locally with 
the Newtonian time ¢ kept by the clock moving with the particle 
considered, it follows that the locally observable properties of the 
Einstein-de Sitter universe are identical with the properties pre- 
dicted for the Newtonian universe. It can be shown that just as in 
the Newtonian case a particle endowed with the parabolic velocity 
steadily decreases in velocity, ultimately to zero, so in the Einstein- 
de Sitter universe the red-shift A’/A, calculated as the Doppler effect, 
for any given particle of the system steadily decreases as the epoch 
of observation of this Doppler effect advances. Thus in the Einstein- 
de Sitter universe, as in the Newtonian universe, each particle may 
be described as undergoing deceleration. This accounts for the 
shorter time-scale as compared with the kinematic theory, where the 
deceleration is zero. 

In practice, ‘local phenomena’, or phenomena ‘close to the ob- 
server’ means phenomena within say 150-200 million light years’ 
distance; they include all phenomena as yet accessible to observa- 
tion. An analyst of Newton’s period who had no data on nebular 
velocities would be unable to estimate the ‘age’ ¢ or present mean 


* De Sitter, loc. cit. p. 180, equation (58), on cancelling c. 
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density p, but he would have been led to predict a non-static universe 
(with either expansion or contraction), to predict a velocity-distance 
proportionality at any one epoch, and to obtain the formula 

p = a*/§ 
connecting density and rate of expansion. Thus he would have 
secured all the results yet capable of observational test. 


aG 


6. On obtaining the above results I communicated them to Dr. 
McCrea. It occurred to both of us, independently, to generalize the 
results so as to give the elliptic and hyperbolic cases, as well as the 
parabolic, on Newtonian mechanics. Actually Dr. McCrea sent me 
his results first. The paper which follows contains features due to 
both of us. 











NEWTONIAN UNIVERSES AND THE CURVATURE 
OF SPACE 


By W. H. McCREA (London) and E. A. MILNE (Ozford) 
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THE present paper investigates the relationship between the universes 
of relativistic cosmology and the universes that can be constructed 
using only Newtonian dynamics, Newtonian gravitation, and New- 
tonian relativity. It is shown that the governing differential equa- 
tions are identical in form in the two cases, and that the locally 
observable phenomena predicted on the two theories are indistin- 
guishable. It is further shown that a space of positive curvature 
corresponds to a Newtonian universe in which every particle has a 
velocity less than the parabolic velocity of escape from the observer, 
a space of negative curvature to one in which every particle has a 
velocity greater than the parabolic velocity. Thus universes of 
positive, zero, and negative curvatures correspond to elliptic, para- 
bolic, and hyperbolic Newtonian universes respectively. These results 
are applied to obtain a number of new results concerning Doppler 
effects and accelerations in relativistic universes. Finally an investi- 
gation due to Lemaitre is briefly discussed. 


1. The results obtained in the foregoing paper can be extended as 
follows to the case in which, on Newtonian mechanics, the velocity 
v is not necessarily equal to the parabolic velocity of escape. Let v, 
the velocity of a particle at distance r from the observer at time ¢, be 
radial in direction, and a function of r and t. The equation of motion, 

Dv GM(r 
Dt = ( ) (1) 
t r 
ov 


may be written — 
. ot 


fon = —{rGpr, (2) 


where p is a function of ¢ only. The equation of continuity may be 
written in the form 

ldp 124 

pdt ror 

1 @ 


— — (rv 
r? Or 


ry) = 0. (3) 


Hence, here, 
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is a function of ¢ only, independent of r. Put 
Idp _ 
p dt 
1 a 


—3F(t); 


then (r2v) = 3F(t), 


which integrates in the form 
ry = PF (t)+-G(d), 
Git) 
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i.e. v= rF(t)+ 


Insert this in (2). Then, since p is a function of ¢ only, 
1 , ?’(t) Git) 2G (t) 
= ‘4 Paci: 
=| rr ++ (rrO+ DFO 

must be a function of ¢ only. This requires G(t) = 0. Hence by 

(5), v= rFit). 

Inserting this in (2) as before, we have 

F'(t)+[ FP = —4xGp. 
Writing (6) in the form 
1dr _ 
rdt 
and integrating it following the motion, we have 
r= fR(t), (8) 
where f defines the particle considered and R(t) is a universal function 
of ¢ satisfying 
148 F(t) = : (9) 
R dt : 
Hence (10) 
where B is a constant. Introducing (9) and (10) in (7), we find 
1 d*R inGB 
= — 


F(t), 


Rd? 


of which the integral is 


a _ §nGB 


nee K 

dt i: 

where K is a constant. Hence, by (9), 
SnGB K 

Fe) = [+a 


= [§aG@p+ Api}, 
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where K = ABi. (14) 
Accordingly by (6), v = r[§aGp+ Api}. (15) 


This is simply the Newtonian integral of motion, for since 


M(r) = §zpr*, 
it may be written 
1.2. GM(r) , A/3\8 
jor = AO) (Ee, (16) 
and of course M(r) is constant following the motion. 

It follows that the particle, and so every particle, possesses the 
elliptic, parabolic, or hyperbolic velocity according as A=0. The 
constant A is the same for all particles. By (6) or (15), v obeys a 
velocity-distance proportionality at any one epoch; and (15) gives 
explicitly the connexion between the mean density p and the coeffi- 
cient in the velocity-distance proportionality. 

Differentiating (16) following the motion, we see that v obeys (1), 
so that (16) is an actual solution. Moreover, it is clear that the solu- 
tion satisfies Einstein’s cosmological principle. The density p is 
obtained as a function of ¢ by integrating (12) in the form 


) 
t } 
=f =. a= | aaanr (17) 
(§7GB+ KR) (§7G+Aé)? 
0 

where p = 1/8. 

2. Comparison with the equations of relativistic cosmology. 
By (10), equation (11) may be written 


2 @R 

Roan — as) 
where Kk = 82rG/c?*. 
Similarly, (12) may be written 


1 dR 
(asa) +3 ; (19) 


where k=—-5=- ; (20) 
c 
Adding (18) and (19) we have 
2 dR 1dR 
— 0. 19’ 
Rede (asa) +O aed 


Equations (19) and (19’) are formally identical with the equations 
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of relativistic cosmology* for an expanding universe of ‘radius’ R, 
curvature k/R?, with A= 0, p= 0. Thus the Newtonian distance 
r is the same function of the Newtonian time ¢ as the ‘distance’ r is 
of ‘cosmic time’ ¢ in the relativistic solutions. Further, we see that an 
expanding space of positive curvature (k > 0) corresponds to a 
Newtonian universe with elliptic velocities (A < 0), an expanding 
space of negative curvature (k < 0) to a Newtonian universe with 
hyperbolic velocities (A > 0).+ By choice of a multiplying factor 
for R, which is equivalent to choice of B, k may be reduced to +1. 
The constant B has no physical significance; it disappears from all 
formulae relating observable quantities. The constant A alone is of 
physical significance.{ 

It follows that the local properties of the universes in expanding 
spaces of positive, zero or negative curvatures are observationally the 
same as in Newtonian universes with velocities respectively less than, 
equal to, or greater than the parabolic velocity of escape. This gives 
great insight into the physical significance of expanding curved space. 


3. The ‘cosmological’ terms involving the cosmical constant A are 
at once obtained, if we superimpose on Newtonian gravitation a 
repulsive force proportional to distance. By a simple application of 
the triangle of forces the equation of motion (1) becomes modified to 

Dv = GM(r) 

Dt }§6orr 
and the equation of continuity is unaffected. The integrations can 
now be carried out as before,§ and we are led to equations for R of 
the relativistic form with the usual A-terms. Einstein’s cosmological 
principle is still satisfied. But from this point of view the introduction 
of A-terms is somewhat artificial. Their introduction is permissible 
simply because they are the only new type of action at a distance 
(in addition to Newton’s law) compatible with the satisfaction of the 
cosmological principle. This last point may be worthy of remark as 
throwing light on the nature of a ‘law of gravitation’. 





jody, (21) 


* e.g. H. P. Robertson, Reviews of Modern Physics, 5 (1933), 62-90, equa- 
tions (3.2). It is clear that (18) and (19) are formally identical with the 
Newtonian equation of motion and its first integral, and that they imply the 
equation of continuity. 

+ Cf. de Sitter, Univ. California Pub. Math. 2 (1933), 171. 

t e.g. the curvature k/R® is simply — Ap! /c’. 

§ The velocity-distance proportionality becomes v = r[$7Gp+ Ap! +}ec*Al}}. 
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4. Omitting further consideration of the A-terms, we may put, in 
(17), 6 = §7@¢/|A|, when it becomes 


t= a f gab 7 
|A| ‘ 

s [i+ r9| (22) 
If A > 0 (hyperbolic case), as ¢ increases, ¢ or @ steadily increases, 
and p steadily decreases to zero. If A < 0 (elliptic case) @ can never 
exceed $7G/|A|, and the density p has a lower limit, after which it 
increases again. We then have an oscillating universe. It should be 
noted that in all cases Dv/Dt is always negative, so that in the 
Newtonian case each particle is steadily decelerated. 

In the hyperbolic case v/r > 0 following any particle, and v > const., 
though to a different constant for each particle corresponding to a 
different value of the parameter f. Also, since p > 0, the asymptotic 
form corresponds to the ‘hydrodynamic’ solution previously obtained* 
by one of us from kinematic principles only. This has been previously 
shown by the methods of general relativity,+ but owing to the formal 
identity of the various relations with the corresponding Newtonian 
ones it holds in the Newtonian case as well. In particular, in this 
asymptotic case we have v = r/t for every particle. We can easily 
show, still of course omitting cosmical repulsion, that in all other 





cases v<rt. (23) 
For the integral (12) may be written 
ee i 

where C = §t@B = ixc?B > 0 since the density is positive. Hence 
for small R we have t~ 80-4RI, (25) 
where we choose ¢ = 0 for R = 0. Also from (6) and (9), 

v= rR'/R. 
Therefore we shall have v < r/t, provided 

t< R/R’. (26) 


* E. A. Milne, Zeits. fiir Astrophys. 6 (1933), 1-95, 71 et seq. 

+ W. O. Kermack and W. H. McCrea, M.N.R.A.S. 93 (1933), 519-29; 
also H. P. Robertson, Zeits. fiir Astrophys. 7 (1933), Heft 3. In the present 
paper we leave aside the question of the wider significance of the kinematic 
solution. 
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But from (24), 
ai = alompEy) = (ors = ate eh 
dt\R’} dt\(CR4+K)]/ (CR-!+K)! CR1+K’ ~ 
and this exceeds unity for K > 0. Thus in the hyperbolic case R/R’ 
increases faster than ¢, and, from (25), R/R’ = 0 when t = 0, so that 
the inequality (26) holds for all ¢ > 0. Thus (23) is established in 
the hyperbolic case. In the elliptic case a similar proof holds so long 
as R’ > 0, and after that the inequality is true a fortiori. This 
shows in a simple manner the retarding effect of gravitation on the 
scattering of the particles. 

It is of some interest to integrate (17), or equivalently (24); (17) is 
more significant, since the meaningless multiplier B is absent. We 
find, putting §7G = a, 


— Oa—A’0)!, a _. _,[A’\3 = 
3 +,sin ( z= (A = —A’ < 0,k > 0), 


t= 


a 


64(a-+A6)! a (= : 


noe 
A Ai 


where it may be recalled that p = 1/4°. 


=) (A>0, k <0), 


a 


5. Doppler effect in curved universes. The general features 
of the motions studied in this paper suggest that in a relativistic 
universe the Doppler effect for any given particle will decrease as 
the epoch of observation advances, and for negative curvatures tend 
to a constant limit, different for different particles. We shall now 
show that this actually is the case. 

We take the metric for the ‘expanding-space’ universe in the form* 

ds? = c? dt?— R? du?, 
so that the particles have fixed coordinates in the space represented 
by du®. Hence if ¢ is the ‘cosmic time’ of departure of a light-signal, 
t, its time of arrival at the observer, we have 


te 


| an =: | du = const., (28) 


* H. P. Robertson, loc. cit. 
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where f du is the coordinate distance along the light-track, for which, 
of course, ds = 0. Hence 
dt, dt 
a ice ae 29 
Rl) RO) = 
If then A is the emitted wave-length, A, the observed wave-length, 
then 
A, _ dt, _ Rite) 
ft ax ates Ge > B, 
dt RO” — 
since R(t) is an increasing function of ¢, and ¢, >¢. This gives the 
red-shift. Consequently, this shift will decrease as the time of obser- 
ration t, adv if 
vation ¢, advances, i d [Rlts) 
dt, \ R(t) 
R'(t, 
2 


)<» 


) R'(t) 
Rie, dts— Fa dt < 0, 


that is, using (29), if R' (tz) < R(t), (31) 
that is, if R’(t) is a decreasing function of t. But we have proved, (cf. 
equation (8)), that R(¢) is the same function of ¢, apart from a con- 
stant multiplier, as the Newtonian distance r is of Newtonian time ¢, 
in the corresponding Newtonian universe. Now in this universe every 
particle is decelerated. Hence dr/dt? < 0, or d*R/dt? <0. Hence 
dR/dt is a decreasing function of t. Therefore the Doppler effect de- 
creases with advancing epoch of observation. Also, in the asymptotic 
form of the hyperbolic case, R(t) = ¢ const., and so 


A,/A = t,/t = const., 


that is, if - 


from (28).* 

These results are worthy of note, since it is sometimes supposed 
that particles necessarily undergo acceleration and not deceleration. 
As a matter of fact acceleration away from the observer is possible 
only when cosmical repulsion predominates over gravitational 
attraction. 


6. Lemaitre’s theory of condensations. A further instance 
of the applicability of the same physical interpretation of the equa- 
tions of general relativity as is given in this paper is provided by 
Lemaitre’s theory of the formation of condensations.t 


* Kermack and McCrea, loc. cit. equation (21). 
+ G. Lemaitre, Comptes Rendus, 196 (1933), 903-4, 1085-7. 
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He studies the motion of a distribution of matter having spherical 
symmetry about the origin, and obtains for any particle an ‘equation 
of motion’ ap\2 _m(x,t) Pe 

(F) = —c*sin*y+ 2G — + $Ac*r?, (32) 
ot r 
where r is its ‘distance’ from the origin, (x, 0, 6) are the coordinates of 
the particle in the space in which it is fixed, and 
r(x, t) 
m(x,t) = 4rpr? dr 
0 

gives a measure of the mass inside radius r. With the addition of the 
cosmical constant A, (32) is the analogy for Lemaitre’s case of our 
equation (16). 

He proceeds to treat the case where m(x,¢) is such a function of 
x alone that there exists a real value x, such that the positive values 
of r for which ér/ét = 0 are imaginary, coincident, or real, according 
as x > xo X = Xo X < Xo It follows that, if r is small when ¢ is 
small, then, if y > xo, r varies from zero to infinity; if x = x», r tends 
asymptotically to a value r, > 0; if x < x, 7 increases up to the first 
positive root of ér/ét = 0, and then decreases again to zero. That is 
to say, the matter outside x = x, continues to spread away from 
the centre, while the matter inside y = y, ultimately falls back upon 
itself. This, then, is the mode of formation of condensations, for 
example, the extra-galactic nebulae, contemplated by Lemaitre. He 
points out that it is necessary for the form of this theory that A + 0. 

Without going further into the implications of Lemaitre’s assump- 
tions, it is now clear that his work would be interpreted in Newtonian 
language by saying that he is dealing with such a distribution of 
matter and motion that all the particles outside a certain shell are 
endowed with a ‘hyperbolic’ radial velocity of escape from the matter 
between themselves and the centre, while all the particles inside this shell 
have an ‘elliptic’ radial velocity, and so must fall back on the centre.* 


* Added in proof: In a very recent paper (Proc. Nat. Acad. Sci. 20 (Jan. 
1934), 12-17), seen after the above was written, Lemaitre has himself alluded 
to the classical analogy to his work, without however following it out in 
detail or showing the connexion with the equation of continuity. 
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